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BACKGROUND OF THE INVENTION 

The present invention relates generally to processes for purification of 
hydrocarbons and, more particularly, to adsorption processes using sorbents selective to 
sulfur compounds and to aromatic compounds. 

Petroleum is an extremely complex mixture and consists predominantly 
of hydrocarbons, as well as compounds containing nitrogen, oxygen, and sulfur. Most 
petroleums also contain minor amounts of nickel and vanadium. The chemical and 
physical properties of petroleum vary considerably because of the variations in 
composition. 

The ultimate analysis (elemental composition) of petroleum tends to 
vary over relatively narrow limits — carbon: 83.0 to 87.0 percent; hydrogen: 10.0 to 
14.0 percent; nitrogen: 0.1 to 1.5 percent; oxygen: 0.1 to 1.5 percent; sulfur: 0.1 to 5.0 
percent; metals (nickel plus vanadium): 10 to 500 ppm. 

Crude oils are seldom used as fuel because they are more valuable when 
refined to petroleum products. Distillation separates the crude oil into fractions 
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equivalent in boiling range to gasoline, kerosene, gas oil, lubricating oil, and residual 
Thermal or catalytic cracking is used to convert kerosene, gas oil, or residual to 
gasoline, lower-boiling fractions, and a residual coke. Petrochemical intermediates 
such as ethylene and propylene are primarily produced by the thermal cracking of light 
5 hydrocarbon feedstocks in the presence of steam. Catalytic reforming, isomerization, 
alkylation, polymerization, hydrogenation, and combinations of these catalytic 
processes are used to upgrade the various refinery intermediates into improved gasoline 
stocks or distillates. The major finished products are usually blends of a number of 
stocks, plus additives. 

10 Gasoline is a complex mixture of hydrocarbons that distills within the 

range 100 to 400 °F. Commercial gasolines are blends of straight-run, cracked, 
reformed, and natural gasolines. Straight-run gasoline is recovered from crude 
petroleum by distillation and contains a large proportion of normal hydrocarbons of the 
paraffin series. Cracked gasoline is manufactured by heating crude-petroleum 

1 5 distillation fractions or residues under pressure, or by heating with or without pressure 
in the presence of a catalyst. Heavier hydrocarbons are broken into smaller molecules, 
some of which distill in the gasoline range. Reformed gasoline is made by passing 
gasoline fractions over catalysts in such a manner that low-octane-number 
hydrocarbons are molecularly rearranged to high-octane-number components. Many of 

20 the catalysts use platinum and other metals deposited on a silica and/or alumina 

support. Natural gasoline is obtained from natural gas by liquefying those constituents 
which boil in the gasoline range either by compression and cooling or by absorption in 
oil. 

Removal of the sulfur-containing compounds is an important operation 
25 in petroleum refining, and is achieved by catalytic processes at elevated temperatures 
and pressures. See, Farrauto, R. J.; Bartholomew, C. H. Fundamentals of Industrial 
Catalytic Processes, Chapman and Hall, New York (1997). The hydrodesulfurization 
(HDS) process is efficient in removing thiols and sulfides, but much less effective for 
heterocyclic diunsaturated sulfur compounds, such as thiophenes and thiophene 
30 compounds/derivatives (e.g., benzothiophene and dibenzothiophene). 
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In 1998, the U.S. largest automakers pledged to put clean-burning cars 
on the road by 2001, beating the Clean Air Act Amendments mandate by five years. 
They proposed the use of internal combustion (IC) engines capable of emitting 50% 
fewer nitrogen oxides (NO x ) and 70% fewer hydrocarbons, thanks to advanced catalytic 
5 converters. Shortly after this low-emission vehicle concept was announced, the U.S. 

Environmental Protection Agency (EPA) revealed concerns that these reductions might 
not be achievable if high-sulfur gasoline and diesel fuel continue to be used. Studies 
involving the EPA and the automobile and oil industries showed that fuel sulfur atoms 
can bond with reactive sites on the catalyst surface, preventing catalyzed reactions 

10 needed to break down NO x and hydrocarbons. Since high-sulfur gasoline may perhaps 
decrease the effectiveness of advanced catalytic converters, the EPA mandates a 
reduction in gasoline and diesel sulfur levels to 30 and 15 ppm, respectively, down 
from the current levels of 300-500 ppm. Krause, C. An Emissions Mission: Solving 
the Sulfur Problem. Oak Ridge National Lab Review (2000), 33, 3. This should be 

1 5 attained by the year 2006. Faced with the severely high costs of compliance, a 
surprising number of refiners are seriously considering reducing or eliminating 
production of on-board fuels. See Parkinson, G., "Diesel Desulfiirization Puts Refiners 
in a Quandary," Chemical Engineering (2001), February issue, 37. 

Ultra-clean fuel may also be desirable for use with a fuel cell system. 

20 For the automotive fuel cells, liquid hydrocarbons may be ideal fuels due to their higher 
energy density, availability, and safety for transportation and storage. However, liquid 
hydrocarbons usually contain certain sulfur compounds that are poisonous to both the 
shift catalysts in the hydrocarbon fuel processors and the electrode catalysts in fuel cell 
processes. Thus, the sulfur content in the liquid hydrocarbons would desirably be 

25 generally less than about 0.1 ppm. 

During the last decade, there have been several published accounts on 
using adsorption for liquid fuel desulfurization. Commercially available sorbents (i.e., 
zeolites, activated carbon and activated alumina) were used in all of these studies. 
Weitkamp et al. reported that thiophene adsorbed more selectively than benzene on 

30 ZSM-5 zeolite. See Weitkamp, J.; Schwark, M.; Ernest, S. "Removal of Thiophene 
Impurities from Benzene by Selective Adsorption in Zeolite ZSM-5," J. Chem. Soc. 



Chem. Commun. (1991), 1 133. Without being bound to any theory, it is believed that 
this is because thiophene (C 4 H 4 S, also known as thiofiiran) has a higher dipole moment 
(0.55 debye) than benzene (non-polar), although their polarizabilities are similar. 
Based on this study, King et al. studied selective adsorption of thiophene, methyl- and 
5 dimethyl-thiophenes (all with one ring) over toluene and p-xylene, also using ZSM-5. 
See King, D. L.; Faz, C; Flynn, T. "Desulfurization of Gasoline Feedstocks for 
Application in Fuel Reforming/' SAE Paper 2000-01-0002, Soc. Automotive 
Engineers, Detroit (2000). They showed that thiophene was more selectively adsorbed, 
both based on fixed bed breakthrough experiments. However, the capacities for 

10 thiophene were unfortunately quite low (only 1-2% wt. adsorbed at 1% thiophene 
concentration). Both vapor phase and liquid phase breakthrough experiments were 
done in these studies, and the results from two phases were consistent. 

The pore dimensions of ZSM-5 are 5.2-5.6 A. Hence, organic sulfur 
compounds with more than one ring will be sterically hindered or excluded. Zeolites 

1 5 with larger pores, as well as larger pore volumes, would appear to be more desirable 
than ZSM-5 as the selective sorbents. Indeed, results of Salem and Hamid indicated 
that 13X zeolite as well as activated carbon had much higher sorption capacities for 
sulfur compounds. See Salem, A. S. H.; Hamid, H. S. "Removal of Sulfur Compounds 
from Naphtha Solutions Using Solid Adsorbents," Chem. Eng. Tech. (1997), 20, 342. 

20 Based on the data of Salem and Hamid, Id., the capacity for sulfur compounds by 13X 
zeolite was approximately an order of magnitude higher than that of ZSM-5, when 
compared with the data of King et al. (cited above) extrapolated to the same conditions. 

Activated alumina (Alcoa Selexsorb) has been used in an adsorption 
process by Irvine. See, for example, U. S. Patent No. 5,730,860 issued to R.L. Irvine in 

25 1998, entitled "Process for Desulfurizing Gasoline and Hydrocarbon Feedstocks." 

No direct comparison has been made among these commercial sorbents. 
Their experiments were mostly done in fixed bed adsorbers, by measuring the 
breakthrough capacities. Based on the literature, the large pore zeolites (NaX or NaY) 
are about the same as activated carbon and alumina, in terms of adsorption of 

30 thiophene. 
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As mentioned hereinabove, current sulfur levels found in commercial 
liquid fuels are commonly obtained by Hydrodesulfiirization (HDS) treatment. As 
further mentioned above, this method is very effective in removing thiols and sulfides, 
but it is generally not adequate for the removal of thiophenic compounds. For instance, 
5 the H 2 S produced during reaction of some thiophene derivatives is one of the main 
inhibitors for deep HDS of unreactive species. Ma, X.; Sakanishi, K.; Mochida, I. 
Hydrodesulfiirization Reactivities of Various Sulfur Compounds in Diesel Fuel. Ind. 
Eng. Chem. Res. (1994), 33, 218 and Knudsen, K.G.; Cooper, B.H.; Topsee, H. 
Catalyst and Process Technologies for Ultra Low Sulfur-Diesel. Appl CataL A-Gen. 

10 (1999), 189, 205. For HDS to meet the new federal government mandates, it is 

believed that reactors with volumes 5-15 times larger (depending on the H 2 pressure) 
than those currently used may be needed. 

Ma et al. studied fixed-bed adsorption of thiophene compounds from jet 
fuels and diesel using an undisclosed transition metal compound (5 wt% loading) 

15 supported in silica gel. Ma, X.; Sun, L.; Song, C. A New Approach to Deep 

Desulfurization of Gasoline, Diesel Fuel and Jet Fuel by Selective Adsorption for Ultra- 
Clean Fuels and for Fuel Cell Applications. Catal. Today (2002), 77, 107 and Ma, X.; 
Sprague, M.; Sun, L.; Song, C. Deep Desulfurization of Liquid Hydrocarbons by 
Selective Adsorption for Fuel Cell Applications. Am. Chem. Soc. Div. Pet. Chem. 

20 Prepr. ( 2002), 47, 48. For jet fuel, they obtained a saturation adsorption capacity of 
0.015 g of sulfur per cm 3 of adsorbent and also showed that breakthrough occurs at 
about 20 cm 3 effluent volume for about 3.2 cm 3 of the metal loaded silica gel. For a 
model diesel fuel, Ma et al. obtained a breakthrough capacity of 1 cm 3 per gram of 
adsorbent. The latter was done for removal of dibenzo thiophene and 4,6- 

25 dibenzothiophene molecules only. Collins et al. also performed fixed-bed adsorption 
experiments for sulfur removal, but after oxidation of the thiophenic compounds. 
Collins, F.M.; Lucy, A.R.; Sharp, C. Oxidative Desulphurisation of Oils via Hydrogen 
Peroxide and Heteropolyanion Catalysis. J. Mol. Catal. A-Chem. (1997), 77 7, 397. 
Oxidation was accomplished by using hydrogen peroxide, an acid catalyst and a phase 

30 transfer agent. Afterwards, the oxidized sulfurs were removed from the diesel oil using 
a silica gel. A breakthrough capacity of about 1 1 cm 3 per gram of silica gel was 
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obtained in this case. Another adsorbent that has been studied was ALCOA Selexsorb, 
which is an activated alumina. In one specific application (Irvine, R.L. Process for 
Desulfurizing Gasoline and Hydrocarbon Feedstocks. U.S. Patent 5,730,860, 1998) this 
proprietary material was used in a temperature swing adsorption (TSA) process in order 
5 to continuously adsorb hetereoatoms from hydrocarbon mixtures and produce full 
boiling range FCC gasoline products with a maximum sulfur content of 30 ppmw. 

Essentially all industrial adsorption processes are based on van der Waals 
interactions between the sorbate and the sorbent. Chemical bonds have yet to be 
exploited. Further, the drawbacks concomitant with HDS appear to render it an 

10 inappropriate solution. 

Thus, it would be desirable to provide an adsorption process for selectively 
removing sulfur compounds from liquid fuels at ambient conditions, thereby 
advantageously leading to a major advance in petroleum refining. It would further be 
desirable to provide highly selective sorbents for this process, thereby overcoming the 

15 drawbacks of current commercial sorbents, which are not desirable for this application. 

SUMMARY OF THE INVENTION 

The present invention addresses and solves the above-mentioned drawbacks, 
and substantially meets the desiderata and advantages enumerated above, by providing 

20 a process for removing thiophene and thiophene compounds from liquid fuel. The 
method comprises the step of contacting the liquid fuel with an adsorbent which 
preferentially adsorbs the thiophene and thiophene compounds, at a selected 
temperature and pressure, thereby producing a non-adsorbed component and a 
thiophene/thiophene compound-rich adsorbed component. The adsorbent comprises 

25 any ion-exchanged zeolite, but in a preferred embodiment, the zeolite is selected from 
the group consisting of zeolite X, zeolite Y, zeolite LSX, and mixtures thereof. The 
zeolite has exchangeable cationic sites, and at least some of the sites have d-block 
transition metal cation present. The preferential adsorption occurs by rc-complexation. 

In a further embodiment, the invention comprises a process for removing 

30 thiophene and thiophene compounds from liquid fuel. The method comprises the step 
of contacting the liquid fuel with an adsorbent which preferentially adsorbs the 
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thiophene and thiophene compounds, at a selected temperature and pressure, thereby 
producing a non-adsorbed component and a thiophene/thiophene compound-rich 
adsorbed component. The adsorbent comprises a carrier having a surface area, the 
carrier having a d-block transition metal compound dispersed on at least some of the 
5 surface area. The metal compound releasably retains the thiophenes. The carrier 

comprises a plurality of pores having a pore size greater than the effective molecular 
diameter of the thiophene/thiophene compounds. The method may further comprise the 
step of changing at least one of the pressure and temperature to thereby release the 
thiophene/thiophene compound-rich component from the adsorbent. 
10 In a further embodiment, the invention comprises a process for selective 

removal of aromatics from hydrocarbons. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Objects, features and advantages of the present invention will become 
15 apparent by reference to the following detailed description and drawings, in which: 

Fig. 1 is a graph depicting equilibrium isotherms of benzene and 
thiophene on NaY (Si/Al = 2.43) at 120° C; with the lines being fittings by Langmuir- 
Freundlich isotherm equation; 

Fig. 2 is a graph depicting equilibrium isotherms on AgY (56 AgVunit 
20 cell) and CuY (14 CuVunit cell) at 120°C; 

Fig. 3 is a graph depicting thiophene/benzene and thiophene/n-octane 
breakthrough in AgY, NaY zeolite at room temperature; 

Fig. 4 is a graph depicting cumulative pore volumes of sorbents used in 

Experiment C; 

25 Fig. 5 is a graph depicting pure component equilibrium isotherms of 

benzene and thiophene on Na-Y (Si/Al=2.43) at 120°C and 180°C; curves are fitted 
with Dubinin- Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; 

Fig. 6 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on Ag-Y (Si/Al=2.43) at 120°C and 180°C; curves are fitted 

30 with Dubinin- Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; 
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fitted curves are not shown for benzene adsorption at 180°C because the artificial 
crossovers to the curves for thiophene at 180° C are observed; 

Fig. 7 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on Cu-Y (Si/Al=2.43) at 90°C and 120°C; curves are fitted with 
5 Dubinin- Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; 

Fig. 8 is a graph depicting pure component equilibrium isotherms of 
thiophene on Ag-X (Si/Al=1.25) and Ag-Y(Si/Al=2.43) at 120°C; 

Fig. 9 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on H-USY (Si/Al=195) at 90°C and 120°C; curves are fitted 
10 with Dubinin- Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; 

Fig. 10 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on Na-ZSM-5 (Si/Al=10) at 90°C and 120°C; curves are fitted 
with Dubinin- Astakhov (solid line) and Langmuir-Freundlich (dotted line) isotherms; 

Fig. 1 1 is a graph depicting pure component equilibrium isotherms of 
15 benzene and thiophene on activated carbon (Type PCB) at 90 °C and 120°C; curves are 
fitted with Dubinin- Astakhov (solid line) and Langmuir-Freundlich (dotted line) 
isotherms; 

Fig. 12 is a graph depicting pure component equilibrium isotherms of 
benzene and thiophene on modified activated alumina (Selexsorb CDX) at 90°C and 
20 120°C; curves are fitted with Dubinin-Astakhov (solid line) and Langmuir-Freundlich 
(dotted line) isotherms; 

Fig. 13 is a graph comparing equilibrium adsorption isotherms of 
thiophene at 120°C; 

Figs. 14-21 are graphs depicting various breakthrough curves; 
25 Fig. 22 is a graph depicting pure component equilibrium isotherms for 

benzene and cyclohexane adsorption on H-USY at 120 and 180°C; curves are fitted 
with Dubinin-Astakhov (solid lines) and Langmuir-Freundlich (dotted lines) isotherms; 

Fig. 23 is a graph depicting pure component equilibrium isotherms for 
benzene and cyclohexane adsorption on Na-Y at 120 and 180°C; curves are fitted with 
30 Dubinin-Astakhov (solid lines) and Langmuir-Freundlich (dotted lines) isotherms; 
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Fig. 24 is a graph depicting pure component equilibrium isotherms for 
benzene and cyclohexane adsorption on Ag-Y at 120 and 180°C; curves are fitted with 
Dubinin-Astakhov (solid lines) and Langmuir-Freundlich (dotted lines) isotherms; 

Fig. 25 is a graph depicting pure component equilibrium isotherms for 
5 benzene and cyclohexane adsorption on Pd-Y at 120°C; curves are fitted with Dubinin- 
Astakhov (solid lines) and Langmuir-Freundlich (dotted lines) isotherms; 

Fig. 26 is a graph depicting pore size distributions of Na-Y, Ag-Y, Pd-Y 

and H-USY; 

Fig. 27 is a graph comparing experimental and GCMC simulation for 
10 adsorption of benzene on Na-Y; 

Fig. 28 is a graph comparing experimental and GCMC simulation for 
adsorption of benzene on H-USY; 

Fig. 29 is a graph comparing experimental and GCMC simulation for 
adsorption of benzene on Ag-Y; 
15 Fig. 30 is a graph depicting thiophene breakthrough curve in reduced 

Cu-Y zeolite at room temperature (C, = 2000 ppmw thiophene; solvent: n-octane); 

Fig. 31 is a graph depicting thiophene breakthrough curve in Ag-Y 
zeolite at room temperature (C 4 - = 2000 ppmw thiophene; solvent: n-octane); 

Fig. 32 is a graph depicting thiophene breakthrough curves in reduced 
20 Cu-Y zeolite at room temperature ((C,= 500 ppmw thiophene; solvent: n-octane); 

Fig. 33 is a graph depicting GC-FPD chromatogram of commercial 
unleaded regular gasoline; 

Fig. 34 is a graph depicting GC-FPD chromatogram of commercial 
diesel gasoline; 

25 Fig. 35 is a graph depicting breakthrough of total sulfur in a fixed-bed 

adsorber with Cu(I)-Y or AC/ Cu(I)-Y adsorbents, during gasoline treatment at room 
temperature. C, is the total sulfur concentration at saturation conditions; 

Fig. 36 is a graph depicting progression GC-FPD chromatograms of 
sulfur compounds during gasoline treatment with Cu(I)-Y. Also shown is processed 

30 fuel volume values normalized by weight of adsorbent; 



10 

Fig. 37 is a graph depicting breakthrough of thiophene or methyl- 
thiophene or benzothiophene or methyl-benzothiophene in a fixed-bed adsorber with 
Cu(I)-Y adsorbent, throughout gasoline treatment at room temperature. C, is the sulfur 
concentration for each species at saturation conditions; 

Fig. 38 is a graph depicting progression GC-FPD chromatograms of 
sulfur compounds during gasoline treatment with AC/ Cu(I)-Y. Also shown is 
processed fuel volume values normalized by weight of adsorbent; 

Fig. 39 is a graph depicting breakthrough of thiophene or methyl- 
thiophene or benzothiophene or methyl-benzothiophene in a fixed-bed adsorber with 
AC/ Cu(I)-Y adsorbent, throughout gasoline treatment at room temperature. C } is the 
sulfur concentration for each species at saturation conditions; 

Fig. 40 is a graph depicting breakthrough of total sulfur in a fixed-bed 
adsorber with AC/ Cu(I)-Y adsorbent, during diesel treatment at room temperature. C, 
is the total sulfur concentration at saturation conditions; 

Fig. 41 is a graph depicting progression GC-FPD chromatograms of 
sulfur compounds during diesel treatment with AC/ Cu(I)-Y. Also shown is processed 
fuel volume values normalized by weight of adsorbent; 

Fig. 42 is a graph depicting breakthrough of dibenzo thiophene (j) or 4- 
methyl-dibenzo thiophene or 4,6-dimethyl-dibenzothiophene or 3,6-dimethyl- 
dibenzothiophene in a fixed-bed adsorber with AC/ Cu(I)-Y adsorbent, throughout 
diesel treatment at room temperature. C, is the sulfur concentration for each species at 
saturation conditions; 

Fig. 43 is a graph depicting FTIR spectra of gasoline or diesel for the C- 
H stretching region; 

Fig. 44 is graph depicting gasoline sulfur breakthrough in CuY zeolite at 
room temperature (flow rate = 0.5 ml/min); and 

Fig. 45 is a graph depicting GC-FPD for desulfurization of gasoline with 

AC/Cu(I)-Y. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention is predicated upon the unexpected and fortuitous 
discovery that thiophene and thiophene compounds are adsorbed slightly more 
selectively via 7i-complexation than is benzene. This is quite counter-intuitive, as it 
5 would be expected that benzene, having more double bonds (3) and more n electrons 
than thiophene (2 double bonds), would be more selectively adsorbed via n- 
complexation. An example of a compound having more double bonds being more 
selectively adsorbed than a compound having fewer double bonds may be found in U.S. 
Patent No. 6,215,037, issued to Padin, Munson and Yang in 2001 and entitled, "Method 

10 for Selective Adsorption of Dienes." 

Without being bound to any theory, it is believed that this counter- 
intuitive, slightly higher selectivity for thiophene/thiophene derivatives/compounds 
may be explained by the following theory. The sulfur atom in the five-member 
thiophene ring has more electrons than the carbons. As such, the S atom, with its 

1 5 available electrons and relatively strong attraction, may be aiding in the rc-complexation 
bonding, thus contributing to the higher selectivity of the present sorbents for 
thiophenes over benzene. How much the sulfur could contribute to 7t-complexation 
bonding, however, is not predictable. 

The present invention is a startling and useful discovery, in that the 

20 sorbents and process of the present invention adsorb more than 1 mmol/gram of 

thiophene at 10" 5 atm. The present inventive process utilizing sorbents which adsorb 
via 7c-complexation may be used to remove sulfur compounds (for example, 
heterocyclic diunsaturated sulfur compounds, including but not limited to thiophene, 
thiophenic compounds, and/or derivatives thereof) from hydrocarbons such as liquid 

25 fuels (including, but not limited to gasoline, diesel fuels, coal and shale derived liquid 
fuels, methanol, and the like). Liquid fuels have a high concentration of aromatics, 
such as, e.g., benzene, and a low concentration of thiophene and thiophene compounds 
(in the parts-per-million (ppm)). Thus, it is one important aspect of the present 
invention that the inventive process adsorbs the thiophenes at least slightly more 

3 0 selectively than benzene. 
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This is a substantial advance over currently available commercial sorbents, 
which adsorb little, if any, thiophene or thiophene compounds at low concentrations. 
For example, activated carbon was recently tested by the present inventors, and was 
found to behave similarly to Na-Y (discussed in further detail below), in that it 
5 adsorbed very little (below 1/10 mmol/gram) benzene and thiophene at 10~ 5 atm. 

Without being bound to any theory, it is believed that the higher sorbent 
capacity of embodiments of the present invention may be due to a sorbent pretreatment 
method of an embodiment of the present invention wherein the sorbent is calcined 
between about 350°C and about 450°C in an inert, dry atmosphere, and then cooled in 

10 a dry atmosphere. The calcining may be carried out for an amount of time ranging 

between about zero hours and about 20 or more hours; or it may be carried out for an 
amount of time ranging between about 5 hours and about 15 hours. The present 
inventors found very satisfactory desulfiirization results using sorbents that had been 
calcined for about 10 hours. 

1 5 The present inventive process may advantageously be run at ambient 

temperature and pressure, which is highly desirable for a variety of reasons. It is much 
less energy consuming to run processes at ambient temperature and pressure. Further, 
conventional catalytic desulfurizing processes, processes are typically run at high 
temperatures, e.g., 700°C, and high pressures, causing the (very expensive) catalyst to 

20 be continually deactivated. Further, the reactor size of a conventional process would 
need to be exponentially increased over present reactor sizes in order to remove 
thiophenes lower than the currently acceptable several hundred ppm. Still further, the 
present inventive process maybe used as a first line desulfurizing process, and/or as a 
clean up desulfurizing process to remove sulfur compounds missed by conventional 

25 processes, such as those processes discussed above (e.g., hydrodesulfurization). It 

could be very useful to have the option to add the present invention as a downstream 
"clean up" desulfurizing process, in that it may not be necessary to revamp current 
refining processes upstream from the present inventive process. This could solve the 
cost-prohibitive problems encountered by some refiners. See Parkinson, G., (cited 

30 above). 
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Without being bound to any theory, it is also believed that the excellent 
desulfurization (of fuels such as, for example, gasoline and diesel) results as discussed 
hereinbelow in embodiments of the present invention may at least in part be due to the 
addition of a relatively thin layer of another sorbent as a guard bed. It is contemplated 
5 that all suitable commercial sorbents may be used as a guard bed. In a non-limitative 
embodiment(s) discussed herein, the present inventors used as the guard bed activated 
carbon, 15% of the bed, added at the inlet of the bed; while the main bed that was doing 
the purification work remained CuY zeolite. This sorbent proved orders of magnitude 
more effective than known sorbent(s). 

10 The sorbents of the present invention, one example of which is Cu(I)Y, 

may be substantially fully regenerated by first air oxidation (e.g., at about 350 °C), 
followed by auto-reduction in an inert atmosphere (e.g., at about 450°C). 

The process and sorbents of the present invention may have quite 
significant cost benefits on the future energy picture. For example, for a full tank of 

15 gas, it is estimated that it would cost only about 70 cents to remove the sulfur from the 
current 350 ppm to below 5 ppm. If some blending were done, to meet the new (2006) 
EPA standard, then it is estimated that it would cost only about 12 cents. These costs 
are based on the assumption that the sorbent is disposable - no regeneration. However, 
it is estimated that the cost could be considerably less if the sorbent were regenerated 

20 (for example, as described above, i.e., air oxidation followed by auto-reduction in inert 
atmosphere). 

Although the process of the present invention has specifically tested Cu-Y 
and Ag-Y (as well as other sorbents), it is to be understood that Type X zeolites may in 
some cases be as good as, or better zeolites than Y zeolites, since more cations are 

25 available in X zeolites. Further, it is to be understood that other zeolites are 

contemplated as being within the scope of the present invention. Still further, it is to be 
understood that all d-block transition metals and/or their ions may be used in place of 
the copper or silver, as it is believed that any d-block element will form 7c-complexation 
bonds with thiophene and thiophene compounds. In particular, it is believed that Ni 2+ 

30 and Pd 2+ would be as effective as Cu + and Ag + . Further, the metals do not need to be 
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ion-exchanged, but rather may be dispersed (monolayer dispersion, island dispersion, 
etc.) on a carrier (such as, for example, silica, alumina, etc.) by any suitable method. 

The present inventors noted that chemical complexation adsorbents, such 
as those for 7r-complexation, have barely been utilized in industrial adsorption 
5 applications. The Tt-complexation bonds are stronger than those formed by van der 
Waals interactions, but they are also weak enough as to be broken by traditional 
engineering means such as increasing temperature and/or decreasing pressure. King, 
C J. Separation Processes Based on Reversible Chemical Complexation. In Handbook 
of Separation Process Technology, Rousseau, R.W., Ed.; Wiley; New York, 1987; 

10 Chapter 15. Therefore, the present inventors have fortuitously tailored and developed 
new adsorbents for processes where selective adsorption is needed, such as in the case 
of sulfur removal from fuels. The present inventors have also fortuitously discovered 
that both cooper and silver exchanged faujasite type zeolites, Ag-Y and Cu(I)-Y, are 
excellent adsorbents for removal of thiophene molecules from liquid hydrocarbon 

15 mixtures. Breakthrough and saturation adsorption capacities obtained for an influent 
concentration of 760 ppmw sulfur (or 2000 ppmw thiophene) in w-octane followed the 
order Cu(I)-Y > Ag-Y > H-Y > Na-Y and Cu(I)-Y > H-Y > Na-Y > Ag-Y, respectively. 
For the case of 190 ppmw sulfur in mixtures containing both benzene and «-octane, 
Cu(I)-Y adsorbed 0.70 and 1.40 wt% sulfur at breakthrough and saturation, 

20 respectively. 

In embodiments of the present invention, 7i-complexation adsorbents are 
used, such as, for example, Cu-Y (auto-reduced) zeolites, to remove sulfur compounds 
from commercial fuels, in particular gasoline and diesel, at room temperature and 
atmospheric pressure using fixed-bed adsorption/breakthrough techniques together with 
25 GC-FPD (Flame Photometric Detector) and FTIR analysis. 

In an embodiment of the present invention, the process for removing 
thiophene and thiophene compounds from liquid fuel comprises the step of contacting 
the liquid fuel with an adsorbent which preferentially adsorbs the thiophene and 
thiophene compounds, at a selected temperature and pressure, thereby producing a non- 
30 adsorbed component and a thiophene/thiophene compound-rich adsorbed component. 
It is to be understood that in the "thiophene/thiophene compound-rich adsorbed 
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component" term as used herein, the thiophenes may not be the greatest amount of 
adsorbate. For example, in the case of liquid fuels having a very high concentration of 
aromatics, and a very low concentration of thiophenes, the adsorbed phase will actually 
contain more benzene than thiophene, because there is much more benzene than 
5 thiophene in the solution. The adsorbed phase may have thiophene/benzene at a ratio 
of about V2, However, that would be enough to advantageously purify the solution. 

The adsorbent comprises any ion-exchanged zeolite, but in a preferred 
embodiment, the zeolite is selected from the group consisting of zeolite X, zeolite Y, 
zeolite LSX, and mixtures thereof. The zeolite has exchangeable cationic sites, and at 

10 least some of the sites have d-block transition metal cation present. The preferential 
adsorption occurs by 7u-complexation. 

In a further embodiment, the invention comprises a process for removing 
thiophene and thiophene compounds from liquid fuel. The method comprises the step 
of contacting the liquid fuel with an adsorbent which preferentially adsorbs the 

1 5 thiophene and thiophene compounds, at a selected temperature and pressure, thereby 
producing a non-adsorbed component and a thiophene/thiophene compound-rich 
adsorbed component. The adsorbent comprises a carrier having a surface area, the 
carrier having a d-block transition metal compound dispersed on at least some of the 
surface area. The metal compound releasably retains the thiophenes. The carrier 

20 comprises a plurality of pores having a pore size greater than the effective molecular 

diameter of the thiophene/thiophene compounds. The method may further comprise the 
step of changing at least one of the pressure and temperature to thereby release the 
thiophene/thiophene compound-rich component from the adsorbent. 

As suggested by King, chemical complexation bonds are generally 

25 stronger than van der Waals interactions (thus giving rise to higher selectivities), yet 

weak enough to be reversible (i.e., to be broken by simple engineering means). See, for 
example, King, C. J., "Separation Processes Based on Reversible Chemical 
Complexation," Handbook of Separation Process Technology, R.W. Rousseau ed., 
Wiley, New York (1987), Chap. 15. Therefore, a tremendous opportunity exists for 

30 developing new sorbents (and new applications in separation and purification) by using 
weak chemical bonds, including various forms of complexation bonds. Based on the 
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principles of 7i-complexation, the inventor(s) have already developed a number of new 
sorbents for a number of applications. These include sorbents for: (a) olefin/paraffin 
separations (see, for example, Yang, R.T.; Kikkinides, E. S. "New Sorbents for Olefin- 
Paraffm Separations by Adsorption via 71-Complexation," AIChE J. (1995), 41, 509; 
5 Rege, S. U.; Padin, J; Yang, R. T. "Olefin-Paraffm Separations by Adsorption: 

Equilibrium Separation by 7t-Complexation vs. Kinetic Separation," AIChE J. (1998), 
44, 799; Huang, H.Y.; Padin, J.; Yang, R. T. "Ab Initio Effective Core Potential Study 
of Olefin/Paraffin Separation by Adsorption via 7c-Complexation: Anion and Cation 
Effects on Selective Olefin Adsorption," J. Phys. Chem. B. (1999), 103, 3206; Huang, 

10 H. Y.; Padin, J.; Yang, R. T. "Comparison of 7r-Complexations of Ethylene and Carbon 
Monoxide and Cu + and Ag + ," Ind. Eng. Chem. Res. (1998), 38, 2720; and Padin, J.; 
Yang, R. T. "New Sorbents for Olefm-Paraffin Separations by Adsorption via 71- 
Complexation: Synthesis and Effects of Substrates," Chem. Eng. Sci. (2000), 55, 
2607.); (b) diene/olefin separation or purification (i.e., removal of trace amounts of 

15 dienes from olefins) (see Padin, J.; Yang, R. T.; Munson, C. L. "New Sorbents for 
Olefin-Paraffm Separations and Olefin Purification for C 4 Hydrocarbons," Ind. Eng. 
Chem. Res. (1999), 38, 3614; Jarayaman, A; Yang, R.T.; Munson, C. L.; Chinn, D. 
"Deactivation of 7c-Complexation Adsorbents by Hydrogen and Rejuvenation by 
Oxidation," Ind. Eng. Chem. Res. (2001); and Takahashi, A.; Yang, R.T.; Munson, 

20 C.L.; Chinn, D. "Cu(I)-Y Zeolite as a Superior Adsorbent for Diene/Olefin Separation," 
Langmuir (2001); and (c) aromatics/aliphatics separation and purification (i.e, removal 
of trace amounts of aromatics from aliphatics) (see Takahashi, A.; Yang, F. H.; Yang, 
R. T. "Aromatics/Aliphatics Separation by Adsorption: New Sorbents for Selective 
Aromatics Adsorption by 7i-Complexation," Ind. Eng. Chem. Res. (2000), 39, 3856; 

25 and Takahashi, A.; Yang, R. T. "New Adsorbents for Purification: Selective Removal 
of Aromatics," AIChE J., (submitted 2001). For purification, vapor phase results can be 
applied directly to liquid phase applications. In the present invention, the inventors 
have used molecular orbital calculations to obtain a basic understanding for the bonding 
between the sorbates and sorbent surfaces, and further, to develop a methodology for 

30 predicting and designing 7i-complexation sorbents for targeted molecules (see, for 

example, Huang, H.Y.; Padin, J.; Yang, R. T. "Ab Initio Effective Core Potential Study 
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of Olefin/Paraffin Separation by Adsorption via 7t-Complexation: Anion and Cation 
Effects on Selective Olefin Adsorption," J. Phys. Chem. B. (1999), 103, 3206; and 
Huang, H. Y.; Padin, J.; Yang, R. T. "Comparison of 7c-Complexations of Ethylene and 
Carbon Monoxide and Cu + and Ag + ," Ind. Eng. Chem. Res. (1998), 38, 2720). 
5 In the present invention, the inventors obtained first results on new 

sorbents for desulfurization using 7t-complexation. CuY and AgY were the sorbents, 
and thiophene and benzene were used as the model system. Since NaX (i.e., 13X) or 
NaY is among the best commercially available sorbents for thiophene, NaY was used as 
the reference for comparison. 

10 The sulfur breakthrough curves with AgY and NaY zeolites (Fig. 3), with 

a liquid feed of thiophene in octane, shows AgY to be an excellent sorbent. 

Further experiments were run with Cu(I)Y. Surprisingly, even though the 
results with AgY were very good, the results with Cu(I)Y were better than those of 
AgY. The sorbent capacity of Cu(I)Y is much higher than that of AgY, and very pure 

15 octane was obtained for extended periods of time. The sulfur capacity of Cu(I)Y zeolite 
was found to be about 2.55 mmol/g (or, 21.4 wt %). This is an extremely high capacity 
because the experiments were removing ppm levels of sulfur. 

Ag-exchanged faujasite was claimed for thiophene removal in U.S. Patent 
No. 4,188,285 issued to Michlmayr. In sharp contrast to the present invention, the '285 

20 preferred temperature for adsorption was 200-350°C, and the sorbent was not 

dehydrated by heat-treatment prior to adsorption. Also in sharp contrast to the present 
invention, the '285 sulfur capacities were very low, at about 0.07-0.15 wt% for Ag-Y. 
In the '285 patent, the highest sulfur capacity (of 0.2 wt%) was obtained with the lowest 
Ag content, with USY zeolite. The '285 sorbent was apparently intended for bonding 

25 the sulfur atom with Ag. It was not intended for 7t-complexation. In contrast, the 

sorbents of the present invention, for example, the Cu(I)Y sorbent has a sulfur capacity 
> 100 times larger. 

Vansant et al. in European Patent Publication No. 0 275 855 investigated 
Cu(II)Y for thiophene removal. In contrast to the present invention, the '855 Cu 2+ 

30 exchanged Y zeolite was purposely heat-treated in air (to 200-550°C) to maintain the 
Cu 2+ in the divalent state. Thus, in contrast to the present invention, in the '855 



publication adsorption by 7i-complexation was not intended. The '855 publication 
discloses low thiophene capacities (the highest was 1.6 wt%). In sharp contrast, the 
thiophene capacities of the present invention were about 21 .4 wt%. 

Figs. 14-21 are graphs depicting various breakthrough curves. 
5 To further illustrate the present invention, the following examples are 

given. It is to be understood that these examples are provided for illustrative purposes 
and are not to be construed as limiting the scope of the present invention. 

EXAMPLES 

10 EXPERIMENT A 

Sorbent Preparation 

Cu + -Y was prepared by ion exchange of Na-Y zeolites (Si/Al=2.43, 56 Al 
atoms/u.c, Strem Chemical) with Cu(N0 3 ) 2 followed by reduction of Cu 2+ to Cu + , since 
Cu(I) is not soluble in water. First, as-received Na-Y was exchanged twice using 

15 excess amounts (10-fold cation-exchange-capacity (CEC) assuming that one Cu 2+ 

compensates two aluminum sites) of 0.5 M Cu(N0 3 ) 2 at room temperature for 24 hours. 
After the exchange, the zeolite suspension was filtered and washed with copious 
amount of de-ionized water. The product was dried at 100°C overnight. Several 
groups have reported reduction of Cu 2+ to Cu + in zeolite in He (i.e., auto-reduction, see 

20 Takahashi, A.; Yang, R.T.; Munson, C.L.; Chinn, D. "Cu(I)-Y Zeolite as a Superior 
Adsorbent for Diene/Olefm Separation," Langmuir (2001); and Larsen, S. C; Aylor, 
A.; Bell, A. T.; Reimer, J. A. "Electron Paramagnetic Resonance Studies of Copper 
Ion-Exchanged ZSM-5," J. Phys. Chem, (1994), 98 y 1 1533). In this study, reduction of 
Cu 2+ to Cu + was carried out in He only at 450 °C for 1 hour. Ag + ion-exchanged Y- 

25 zeolite (Ag-Y) was prepared at room temperature for 24 hours in the same manner as 
Cu 2+ exchange, using 5-fold excess AgN0 3 (0.1M). 

Single component isotherms for benzene and thiophene were measured at 
120°C using standard gravimetric methods. A Shimadzu TGA-50 automatic recording 
microbalance was employed. Helium (Pre-purified grade, Metro welding, 99.995%) 

30 was used as the carrier gas and was first passed through two consecutive gas-wash 

bottles (to ensure saturation), which contained benzene (HPLC grade, Aldrich, 99.9+%) 
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or thiophene (Aldrich, 99+%). After diluting the concentration to the desired value by 
blending with additional helium, the mixture was directed into the microbalance. 



Molecular Orbital Computational Details 

5 Molecular orbital (MO) studies on the 7i>complexation bonding for 

benzene and sorbent surfaces had been investigated recently. See Takahashi, A.; Yang, 

F. H.; Yang, R. T. "Aromatics/Aliphatics Separation by Adsorption: New Sorbents for 
Selective Aromatics Adsorption by 7t-Complexation," Ind. Eng. Chem. Res. (2000), 39, 
3856. In this work, similar MO studies were extended to thiophene and sorbent 

10 surfaces. The Gaussian 94 Program in Cerius2 molecular modeling software (see 
Bowie, J. E., Data Visualization in Molecular Science: Tools for Insight and 
Innovations Addison- Wesley Pub. Co.: Reading, Mass., (1995); chapter 9) from 
Molecular Simulation, Inc. was used for all calculations. See Frisch, M. J.; Trucks, 

G. W.; Schlegel, H. B.; Gill, P. M.W.; Johnson, B. G.; Robs, M. A.; Cheeseman, J. R.; 
15 Keith, T.; Petersson, G. A.; Montgomery, J. A.; Raghavachari, K.; Al-Laham, M. A.; 

Zakrzewski, V. G.; Ortiz, J. V.; Foresman, J. B.; Peng, C. Y.; Ayala, P. Y.; Chen, W.; 

Wong, M. W.; Andres, J. L.; Replogle, E. S.; Gomperts, R,; Martin, R. L.; Fox, D. J.; 

Binkley, J. S.; DeFrees, D. J.; Baker, J.; Stewart, J. P.; Head-Gordon, M.; Gonzalez, C; 

Pople, J. A. Gaussian 94, Revision B.3, Gaussian, Inc., Pittsburgh, PA. (1995). MO 
20 calculations for thiophene and sorbent surfaces were performed at the Hartree-Fock 

(HF) and density functional theory (DFT) level using effective core potentials (ECPs). 

See Hay, P. J.; Wadt, W.R. "Ab Initio Effective Core Potential for Molecular 

Calculations. Potentials for the Transition Metals Atoms Sc to Hg, M J. Chem. Phys. 

(1985), 82, 270; Wadt, W. R.; Hay, P. J. "Ab Initio Effective Core Potential for 
25 Molecular Calculations: Potentials for Main Group Elements Na to Bi," J. Chem. Phys. 

(1985), 92, 284; and Hay, P. J.; Wadt, W. R. "Ab Initio Effective Core Potential for 

Molecular Calculations: Potentials for K to Au Including the outermost core orbitals," 
Chem. Phys. (1985), 82, 299. 

The LanL2DZ basis set (see Russo, T. V.; Martin, R. L.; Hay, P. J. 
30 "Effective Core Potentials for DFT Calculations," J. Phys. Chem. (1995), 99, 17085) is 

a double-^ basis set containing ECP representations of electrons near the nuclei for 
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post-third-row atoms. The reliability of this basis set has been confirmed by the. 
accuracy of calculation results as compared with experimental data. Therefore, the 
LanL2DZ basis set was employed for both geometry optimization and natural bond 
orbital (NBO) analysis. 

5 The restricted Hartree-Fock (RHF) theory at the LanL2DZ level basis set 

was used to determine the geometries and the bonding energies of thiophene on AgCl 
and CuCl. The simplest models with only a single metal chloride interacting with a 
thiophene molecule were chosen for 7t-complexation studies. The optimized structures 
were then used for bond energy calculations according to the following expression: 

10 

^ads — -^adsorbate ^adsorbent ^adsorbent-adsorbate (0 

where E adsorbate is the total energy of thiophene, E adsorbent is the total energy of the bare 
adsorbent, i.e., the metal chloride, and E adsorbent . adsorbate is the total energy of the 
1 5 adsorbate/adsorbent system. A higher value of E ads corresponds to a stronger 
adsorption. 

Natural Bond Orbital (NBO) 

The optimized structures were also used for NBO analysis at the 
20 B3LYP/LanL2DZ level. The B3LYP (see Becke, A. D. "Density Functional 

Thermochemistry, m. The Role of Exact Exchanges," J. Chem. Phys. (1993), 98, 
5648) approach is one of the most useful self-consistent hybrid (SCH) approaches (see 
Becke, A. D.; "A New Mixing of Hartree-Fock and Local Density-Functional 
Theories," J. Chem. Phys. (1993), 98, 1372), it is Beck's 3 -parameter nonlocal exchange 
25 functional (see Becke, A. D.; "Density-Functional Thermochemistry. EL The Effect of 
the Perdew-Wang generalized-gradient correlation correction," J. Chem. Phys. (1992), 
97, 9173) with nonlocal correlation functional of Lee, Yang and Parr. See Lee, C; 
Yang, W.; Parr, R.G. "Development of the Colle-Salvetti Correlation-Energy Formula 
into a Functional of the Electron Density," Phys. Rev. (1988), B37, 785. 
30 The NBO analysis performs population analysis that pertains to localized 

wave-function properties. It gives a better description of the electron distribution in 
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compounds of high ionic character, such as those containing metal atoms. See Reed, A. 
E.; Weinstock, R.B.; Weinhold, F. "Natural Population Analysis," J. Chem. Phys. 
(1985), 83 (2), 735. It is known to be sensitive for calculating localized weak 
interactions, such as charge transfer, hydrogen bonding and weak chemisorption. 
5 Therefore, the NBO program (Glendening, E. D.; Reed, A. E.; Carpenter, J. E., 

Weinhold, F. NBO Version 3.1 (1995)) was used for studying the electron density 
distribution of the adsorption system. 

Results and Discussion 

For purification, as the inventors have shown for the removal of ppm 
levels of dienes from olefins (Takahashi, A.; Yang, F. H.; Yang, R. T. 
"Aromatics/Aliphatics Separation by Adsorption: New Sorbents for Selective 
Aromatics Adsorption by 7i-Complexation," Ind. Eng. Chem. Res., (2000), 39, 3856; 
and Takahashi, A.; Yang, R. T. "New Adsorbents for Purification: Selective Removal 
of Aromatics," AIChEJ., (submitted 2001), one parameter is the separation factor; and 
the separation factor should be evaluated at low concentrations of the solute but at high 
concentrations of the solvent. In other words, the capacity for the dilute solute (at ppm 
level) is an important factor in determining the purification. The separation factor (a) is 
given by: 

av-XMOWi) (2) 

where X and Y are mole fractions in the adsorbed and fluid phases, respectively. In our 
case, component 1 is thiophene and component 2 is benzene. Hence, the ability of the 
25 sorbent for purification depends on the amount of thiophene adsorbed at the ppm levels. 

Figure 1 shows the equilibrium isotherms of benzene and thiophene on 
NaY at 120°C. NaY does not have a selectivity for thiophene; however, it adsorbs both 
benzene and thiophene quite strongly, as evidenced by the still measurable amounts 
adsorbed at partial pressures as low as 10-100 ppm of 1 atm. As expected, a fixed-bed 
30 adsorber experiment using NaY would yield a breakthrough curve (Yang, R. T. Gas 
Separation by Adsorption Processes, Butterworth, Boston (1987), Imperial College 



10 



15 



20 
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Press, London (1997)), as seen in the literature (Salem, A. S. H.; Hamid, H. S. 
"Removal of Sulfur Compounds from Naphtha Solutions Using Solid Adsorbents," 
Chem. Eng. Tech. (1997), 20, 342). The data were fitted by the Langmuir-Freundlich 
isotherm (Takahashi, A.; Yang, R.T.; Munson, C.L.; Chinn, D. "Cu(I)-Y Zeolite as a 
5 Superior Adsorbent for Diene/Olefin Separation," Langmuir (2001); Takahashi, A.; 
Yang, F. H.; Yang, R. T. "Aromatics/Aliphatics Separation by Adsorption: New 
Sorbents for Selective Aromatics Adsorption by 7i-Complexation," Ind. Eng. Chem. 
Res. (2000), 39, 3856; and Takahashi, A.; Yang, R. T. "New Adsorbents for 
Purification: Selective Removal of Aromatics," AIChE J. f (submitted 2001), as shown 

10 by the lines in the figure. 

Figure 2 shows the isotherms of benzene and thiophene on AgY and CuY. 
Note that: (a) compared with NaY, AgY and CuY adsorb both thiophene and benzene 
more strongly at below 10 3 atm, due to 7i-complexation; and (b) on a per-cation basis, 
Cu + shows stronger interactions than Ag + . Compared with Figure 1, these sorbents 

15 adsorbed significantly more thiophene and benzene than NaY at pressures below 10" 3 
atm, and nearly the same amounts at high partial pressures. This result was a clear 
indication of 7t-complexation with Ag + and Cu + ; Na + could not form 7i-complexation 
bonds. From the data in Figures 1 and 2, the separation factors were calculated using 
mixed Langmuir-Freundlich isotherms (Taka), and are shown in Table 1 . However, 

20 these separation factors did not reflect the relative strengths of 7t-complexation between 
Cu + and Ag + because the Cu + exchange was not complete. Nonetheless, the large 
separation factors indicate that AgY and CuY are excellent sorbents for purification 
purposes. 

25 Table 1. Separation factors of thiophene over benzene estimated from pure-component isotherms 
by extended Langmuir-Freundlich Equation at the pressure of 3x1 0" 5 atm for thiophene and 1 atm 
for benzene. 



Adsorbent 


Ag-Y 


Cu-Y 


Na-Y 


Separation Factor 


2800 


50 


0.10 



The neutron activation analyses of the sorbent samples showed that the 
Ag + exchange was 100% but the Cu + exchange was only 23%. On a per-cation basis, 
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the 7t-complexation with Cu + was stronger than that with Ag + . This was indeed 
confirmed by our molecular orbital calculations, as will be discussed. To understand 
the relative strengths of 7C-complexation between Ag + and Cu + , the thiophene adsorption 
amounts at 2 x 10" 5 atm were normalized by Ag + or Cu + content, and the results are 
5 shown in Table 2. It is seen that Cu + could adsorb higher thiophene adsorption amounts 
per cation. In fact, 0.92 thiophene molecule per Cu+ was obtained at 2x10-5 atm at 
120°C. This amount was due to Cu + since the amount adsorbed by NaY at the same 
pressure was negligible. At the same pressure, only 0.42 thiophene/Ag + was obtained. 
This result indicated strong rc-complexation bonds between both Cu + and Ag + , and that 
10 the bond with Cu + was stronger. Further work is in progress to increase the Cu + 

contents on various X/Y and other zeolites and the results will be reported shortly. 



Table 2. Thiophene adsorption amount normalized by cation content (thiophene/cation). 



Adsorbent 


Ag-Y 


Cu-Y 


Cation for n-complexation 


Ag + 


Cu + 


Cation Content in Zeolite (wt%) 


38.2 


8.29 


Amount Adsorbed (thiophene 
molecule/cation) 


0.42 


0.92 



20 The heats of adsorption for 7i-complexation were estimated by taking the 

values at low loadings, and they were obtained from the temperature dependence of the 
isotherms at low pressures. These values are given in Table 3. 



Table 3. Heats of adsorption (kcal/mol) for n-complexation. 



Zeolite 


Thiophene 


Benzene 


Cu(I)-Y 


22.3 


21.8 


Ag-Y 


21.4 


20.1 



30 



Bond Energies, Geometries and NBO results 

The energies of adsorption calculated using equation 1 for thiophene and 
benzene are summarized in Table 4. 
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Table 4. Summary of energies of adsorption for thiophene and benzene in kcal/mol. 



MCI 


E, ds (Thiophene) 


E. ds (Benzene) 


CuCl 


13.5 


12.4 


AgCl 


9.0 


8.6 I 



5 

The theoretical calculations indicate that the 7c-complexation strengths follow the order 
CuCl > AgCl and more importantly, thiophene > benzene. This trend is in agreement 
with the experimental data, in Table 3. Chloride was used as the anion in the 
theoretical calculations, while zeolite framework was the anion in the experiment. It is 

1 0 known from our previous work that the anion has a large effect on the 7c-complexation 
bonds. See, for example, Huang, H.Y.; Padin, J.; Yang, R. T. "Ab Initio Effective Core 
Potential Study of Olefin/Paraffin Separation by Adsorption via 7i-Complexation: 
Anion and Cation Effects on Selective Olefin Adsorption," J. Phys. Chem. B. (1999), 
103, 3206; and Huang, H. Y.; Padin, J.; Yang, R. T. "Comparison of 7i-Complexations 

15 of Ethylene and Carbon Monoxide and Cu + and AgV Ind. Eng. Chem. Res. (1998), 38, 
2720. In the optimized structures of thiophene-MCl complexes, the distance between 
J the thiophene molecule and Cu ion is about 0.3 A shorter than that of thiophene and Ag 
ion. The NBO analysis is summarized in Tables 5 and 6. There is some donation of 
electron charges from the rc-orbital of thiophene to the vacant s orbital of metals known 

20 as o donation and, simultaneously, back-donation of electron charges from the d 

orbitals of metals to 7i*orbital of thiophene or n back-donation. The n back-donation 
appears to be important in all cases. The charge transfer results again confirmed the 
experimental data that the relative strengths of the 7i-complexation bonds follow the 
order: thiophene > benzene and Cu + > Ag + . 

25 



Table 5. Summary of NBO analysis* of /r-complexation between thiophene and MCI. 



MCI 


C-*M interaction 
( a donation) 


M-*C interaction 
(d - 7t * back-donation) 

q 2 


Net Change 

qi + q 2 


CuCl 


0.037 


-0.022 


0.015 


AgCl 


0.022 


-0.014 


0.008 



q x is the amount of electron population change in valence s orbitals of the metal, and ^ is the total L 
amount of electron population decrease in valence d orbitals of the metal. 
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Table 6. Summary of NBO analysis* of /r-complexation between benzene and MCL 



MCI 


C-+M interaction 
(o donation) 

qi 


M-*C interaction 
(d - 7T * back-donation) 


Net Change 

qi + q 2 


CuCl 


0.011 


-0.013 


0.002 


AgCl 


0.003 


-0.007 


0.004 



q t is the amount of electron population change in valence s orbitals of the metal, and is the total 



amount of electron population decrease in valence d orbitals of the metal. 
EXPERIMENT B 

Thiophene Removal Using Sodium and Silver Cation in Faujasite Zeolites 

Work has been done on different adsorbents to remove sulfur from 
mixtures, such as, for example, thiophene from benzene or n-octane based solutions. 
These adsorbents include sodium and silver forms of faujasite type zeolites, 
particularly, Y type zeolites. Initial studies were done for the case of an inlet 
concentration of 2000 ppmw of thiophene in either benzene or n-octane. 

Sorbent Preparation 

Sodium Y Zeolite powder was obtained from Strem Chemicals. Silver Y 
Zeolite was obtained by ion exchanging NaY at room temperature. The solution was 
prepared using a silver nitrate salt and deionized water (-0.2 M). A 4-fold excess 
amount of silver was used to ensure exchange completion. After ion exchanging for 48 
hours, the adsorbent was recovered by filtration and washed with large amounts of 
deionized water. Drying was performed at room temperature. 

Breakthrough Experiments 

Breakthrough experiments were performed in a custom made quartz 
reactor equipped with a glass frit. The adsorbent was loaded inside the reactor (usually 
1 or 2 grams) and activated in situ under a helium atmosphere at 350 °C for 24 hrs. The 
adsorbent was then cooled down to room temperature. At this point, the sorbent was 
wetted by pouring down sulfur free solvent (benzene or n-octane) for about 30 minutes. 
This stream was then changed to one containing 2000 ppmw thiophene. Samples were 
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collected at different time intervals, and the thiophene outlet concentration was 
analyzed using a GC unit equipped with a polar column at 65 °C. 



Preliminary Results 

5 Figure 3 shows breakthrough curves for thiophene in either benzene or 

n-octane in sodium or silver forms of Y Zeolite. 

Thiophene adsorption amounts after saturation. 



NaY AgY 
10 

Saturation Amount (mmol/g) 

Solvent 

Benzene 0.102 0.171 

n-Octane — 0.898 

15 



Integration of the area above the curves yield the adsorption saturation 
amounts presented in the table shown above. For the case of thiophene/benzene 
mixtures, the calculated thiophene adsorption amounts are not that high. Very recent 

20 single component vapor phase adsorption experiments have shown that NaY and AgY 
may adsorb benzene and thiophene in similar quantities. Thus, benzene may compete 
for many of the adsorption sites, perhaps leaving thiophene with fewer available sites. 
However, despite the possibility of competitive adsorption, the presence of silver 
cations results in an increase in some adsorption capacity. For the case of thiophene/n- 

25 Octane mixtures, AgY Zeolite shows great selectivity towards thiophene. Without 

being bound to any theory, it is believed that this indicates that the interaction between 
the silver cations and the thiophene molecules is higher than that of n-Octane 
molecules. 



30 EXPERIMENT C 

In this experiment, the known commercial sorbents such as Na-Y, Na- 
ZSM5, H-USY, activated carbon and activated alumina (Alcoa Selexsorb) were 
included, and a direct comparison was made with Cu-Y and Ag-Y which were the 



sorbents with 7i-complexation capability, 
model system for desulfurization. 
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Thiophene and benzene were used as the 



Sorbent preparation 

5 Various kinds of sorbents were investigated in this work. Four as-received 

sorbents: Na-type Y-zeolite (Na-Y, Si/Al=2.43, 56 Al atoms/unit cell, Strem Chemical), 
H-type ultra-stable Y-zeolite (H-USY, Si/Al=195, 0.98 Al atoms/u.c, HSZ-390HUA, 
TOSOH Corporation), activated carbon (Type PCB, Calgon Carbon Corporation) and 
modified activated alumina (Selexsorb CDX, Alcoa Industrial Chemical), were used in 

10 this study. According to the product datasheets, Selexsorb CDX is formulated for 
adsorption of sulfur-based molecules, nitrogen-based molecules, and oxygenated 
hydrocarbon molecules. Na-Y and H-USY were in powder form (binderless). Since 
activated carbon was in granular form and activated alumina was in pellet form, they 
were crushed into powder form for evaluation. 

1 5 Cu + -Y and Ag-Y were prepared as in Experiment A, above. 1 3X 

(Si/Al=1.25, Linde, lot#945084060002) was used for the preparation of Cu-X (10 fold 
CEC solution of Cu(N0 3 ) 2 , ion-exchanged at 65 °C for 24 hrs, three times) and Ag-X 
(5-fold CEC solution of AgN0 3 , ion-exchanged at RT for 24 hrs, twice). Na-type 
ZSM-5 (Na-ZSM-5) was prepared at room temperature by Na + -exchange of NH 4 -ZSM- 

20 5 (Si/Al=10, SM-24, ALSI-PENTA Zeolite GmbH). 

Isotherm and uptake rate measurements 

The objective of this study is to compare the strength of adsorptive 
interaction between adsorbents and thiophene/benzene. Extremely low partial pressures 

25 at less than 10 5 atm would be necessary to meet this objective if isotherms were 

measured at ambient temperature, because the isotherms at ambient temperature are 
fairly flat and are difficult to compare each other. However, it is very difficult to obtain 
and control such low partial pressures experimentally. Therefore, single component 
isotherms for benzene and thiophene were measured at 90°C, 120°C and 180°C using 

30 standard gravimetric methods. A Shimadzu TGA-50 automatic recording microbalance 
was employed. Helium (Pre-purified grade, Metro welding, 99.995%) was used as the 
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carrier gas and was first passed through two consecutive gas-wash bottles (to ensure 
saturation), which contained benzene (HPLC grade, Aldrich, 99.9+%) or thiophene 
(Aldrich, 99+%). After diluting the concentration to the desired value by blending with 
additional helium, the mixture was directed into the microbalance. 
5 Isosteric heats of adsorption were calculated using the Clausius-Clapeyron 

equation from isotherms at different temperatures. Nitrogen isotherms at 77K 
measured with a Micromeritics ASAP 2010 system were used for pore size distribution 
and pore volume determination. Pore size distributions were calculated with Horvath- 
Kawazoe equation. See Horvath, G., and K. Kawazoe, "Method for Calculation of 

10 Effective Pore Size Distribution in Molecular Sieve Carbon," J. Chem. Eng. Japan 
(1983), 16, 470; Saito, A, and H. C. Foley, "Curvature and Parametric Sensitivity in 
Models for Adsorption in Micropore," AIChE J (1991), 37, 429; and Cheng, L. S., and 
R.T. Yang, "Improved Horvath-Kawazoe Equations Including Spherical Pore Models 
for Calculating Micropore Size Distribution," Chem. Eng. Sci. (1994), 49, 2599. Slit 

1 5 pore model was used for activated carbon, and cylindrical pore model was used for 
ZSM-5 and activated alumina, while spherical pore model was applied to Y-zeolites 
and X-zeolites. As for the parameters for H-K equation such as polarizability, magnetic 
susceptibility, density, etc., default values in Micromeritics ASAP system were used for 
the calculation of pore sizes. See "Accelerated Surface Area and Porosimetry System 

20 Operator's Manual V3.02," Micromeritics Instrument Corp. (1997). The diffusion time 
constants, D/r 2 ( s" 1 ), were calculated from the uptake rates. See Yeh, Y. T, "Diffusion 
and Adsorption of Gases in Molecular Sieves," Ph.D. Dissertation, University of New 
York at Buffalo, Buffalo, New York, (1989). In this experiment, short time region (up 
to 30 % uptake) and spherical adsorbent model were used. 

25 

Chemical Analysis 

The compositions of Cu-Y and Ag-Y were characterized using neutron 
activation analysis (NAA) in the research nuclear reactor of the Phoenix Memorial 
Laboratory at the University of Michigan. The sample was irradiated sequentially for 
30 one minute at a core-face location with an average thermal neutron flux of 2 x 10 12 

n/cm 2 /s. Two separate gamma-ray spectra were then collected for each sample with a 
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high resolution germanium detector: one after 13 -minute decay to determine the 
concentrations of Al, Ag and Cu, and a second after a 1 hour and 56-minute decay to 
analyze for Na. Gamma energy lines at 1779 keV, 632.99 keV, 1039.20 keV and 
1368.6 keV were used for the determination of Al, Ag, Cu and Na concentration, 
5 respectively. 

Molecular Orbital Computational Details and Natural Bond Orbitals are 
as described in Experiment A above. 

Models for Ag-Zeolite (AgZ) and Cu-Zeolite (CuZ) 

10 The zeolite models selected for this study are similar to the ones used by 

Chen and Yang (see Chen, N. and R.T. Yang, "Ab Initio Molecular Orbital Study of 
Adsorption of Oxygen, Nitrogen, and Ethylene on Silver-Zeolite and Silver Halides," 
Ind. Eng. Chem. Res. (1996), 35, 4020), with the molecular formula of (HO) 3 Si- O - 
Al(OH) 3 , and the cation Ag+ or Cu + sits 2 - 3 A above the bridging oxygen between 

15 Si and Al. This is a good cluster model representing the chemistry of a univalent cation 
bonded on site II (SII) of the faujasite framework (Z). Once the optimized structures of 
AgZ and CuZ are obtained at the B3LYP/LanL2DZ level, then a molecule of thiophene 
(C 4 H 4 S) or benzene (C 6 H 6 ) is added onto the cation of the zeolite model, and the 
resulting structure is further optimized at the B3LYP/LanL2DZ level. 

20 Results and Discussion 

Characterization of sorbents. Neutron activation analysis showed that 
Ag/Al and Na/Al ratios in Ag-Y were 1.13 and 0.01, respectively. This is because Ag 4 " 
is known to have higher selectivity to cation sites in zeolites compared to Na + . See 
Breck, D.W. Zeolite Molecular Sieves: Structure, Chemistry and Use; Wiley: New 

25 York (1974). More than 100% Ag ion-exchange ratio was obtained in Ag-Y, because 
some Ag ions were located outside the charge-compensating sites. On the other hand, 
the Al, Cu, and Na contents in Cu-Y were 6.10 ± 0.25 wt%, 6.65 ± 0.05 wt%, and 1.50 
±0.12 wt%, respectively. This means that the Cu/Al molar ratio was 0.463, and the 
Na/Al ratio was 0.289, indicating less selectivity than Ag + . Therefore, Cu ion- 
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exchange for X-zeolite was carried out three times at 65 °C (instead of twice at RT for 
Cu-Y) for 24 hrs to enhance the ion-exchange. 

Figure 4 shows the comparison of the cumulative pore volumes of the 
sorbents used. Activated carbon has the largest pore volume with small pores, while 
5 Na-ZSM-5 has the smallest pore volume. Selexsorb CDX also shows small pores. 

Because silver is heavier than copper and the silver content in Ag-Y (Ag/Al=1.13) was 
more than copper content in Cu-Y (Cu/Al = 0.463), the density of Ag-Y was higher by 
36% (Na 0 . 084 Cu 0J35 Al 0>292 Si 0 . 708 O 2 ; MW = 70.27 vs. Ag 0 .3 30 Al 0 . 292 Si 0 . 708 O 2 ; MW = 95.36). 
Consequently, the pore volume of Ag-Y was about 25 % smaller than that of Cu-Y. 

10 Benzene/Thiophene adsorption isotherms: Na-Y. Figure 5 shows the 

equilibrium isotherms of benzene and thiophene on Na-Y at 120°C and 180°C. The 
data were fitted by the Langmuir-Freundlich and Dubinin- Astakhov isotherms (see 
Padin, J.; Yang, R. T.; Munson, C. L., "New Sorbents for Olefin-Paraffin Separations 
and Olefin Purification for C 4 Hydrocarbons," Ind. Eng. Chem. Res. (1999), 38, 3614; 

15 Jayaraman, A; Yang, R. T.; Munson, C. L.; Chinn, D., "Deactivation of 7t- 

Complexation Adsorbents by Hydrogen and Rejuvenation by Oxidation," Ind. Eng. 
Chem. Res. (2001), 40, 4370; Takahashi, A.; Yang, R. T.; Munson, C. L.; Chinn, D., 
"Cu(I)-Y Zeolite as a Superior Adsorbent for Diene/Olefin Separation," Langmuir 
(2001), 17, 8405; Takahashi, A.; Yang, F. H.; Yang, R. T., "Aromatics/Aliphatics 

20 Separation by Adsorption: New Sorbents for Selective Aromatics Adsorption by 7i- 

Complexation," Ind Eng. Chem. Res. (2000), 39, 3856; and Takahashi, A.; Yang, R. 
T., <e New Adsorbents for Purification: Selective Removal of Aromatics," AIChE J. t 
(submitted 2001), which are shown by the lines in Figure 5. Barthomeuf and Ha 
measured benzene adsorption on NaH-Y (Si/Al=2.43, 81%Na) at 172°C. See, 

25 Barthomeuf, D. and Ha, B.-K, "Adsorption of Benzene and Cyclohexane on Faujasite- 
type Zeolites, Part-1, Thermodynamics Properties at Low Coverage,"/. Chem. Soc, 
Faraday Trans. 1 (1973), 69, 2147. Their isotherm and calculated heat of adsorption 
agreed with the data shown in Figure 5 very well, substantially ensuring the accuracy of 
the isotherms in this work. More benzene was adsorbed on Na-Y than thiophene at 

30 pressures lower than 10' 2 atm. This result indicated that the contribution of the higher 
polarizability of benzene (10.3xl0 24 cm 3 for benzene vs. 9.7xl0" 24 cm 3 for thiophene) 
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was larger than the contribution of the dipole moment (0 debye in benzene vs. 0.55 
debye in thiophene) and lower magnetic susceptibility (9.1x10 29 cm 3 /molecules for 
benzene vs. 9.5x1 0" 29 cm 3 /molecules for thiophene). At higher pressures over 1 x 10 2 
atm, the benzene adsorption amount became smaller owing to the larger molar volume 
5 of benzene than thiophene. Na-Y does not have a selectivity for thiophene; however, it 
adsorbs both benzene and thiophene quite strongly, as evidenced by the still measurable 
amounts adsorbed at partial pressures as low as 10" 4 to 10~ 3 atm. 

Ag-Y and Cu-Y. Figure 6 and Figure 7 show the isotherms of benzene 
and thiophene on Ag-Y and Cu-Y. Compared with Figure 5, these sorbents adsorbed 

10 significantly more thiophene/benzene than Na-Y at pressures below 10" 3 atm, and nearly 
the same amounts at high partial pressures. This result was a clear indication of n- 
complexation with Ag + and Cu + ; since Na + could not form rc-complexation bonds. 
However, the difference in the amounts of thiophene/benzene adsorbed did not reflect 
the relative strengths of 7i-complexation between Cu + and Ag + because the Cu + 

15 exchange was not complete. The neutron activation analyses of the sorbent samples 
showed that the Ag + exchange was 100% but the Cu + exchange was only 46%. 
According to the EPR analysis (described in Takahashi, A., Yang, R. T., Munson, C. L. 
and Chinn, D., "Cu(I)-Y Zeolite as a Superior Adsorbent for Diene/Olefm Separation," 
Langmuir (2001), 17, 8405), only about a half of the Cu 2+ was auto-reduced to Cu + after 

20 our heat treatment at 450°C for lhr in He. On a per-cation basis, the 7t-complexation 

with Cu + was stronger than that with Ag + . This was indeed confirmed by our molecular 
orbital calculations, as will be discussed. To understand the relative strengths of n- 
complexation between Ag + and Cu + , the thiophene adsorption amounts at 2 x 10" 5 atm 
were normalized by Ag + or Cu + content, and the results are shown in Table 7. It is seen 

25 that Cu + could adsorb higher thiophene adsorption amounts per cation. In fact, 0.92 
thiophene molecule per Cu + was obtained at 2xl0" 5 atm at 120°C. This amount was 
due to Cu + since the amount adsorbed by NaY at the same pressure was negligible. At 
the same pressure, only 0.42 thiophene/ Ag + was obtained. This result indicated strong 
7c-complexation bonds between both Cu + and Ag + , and that the bond with Cu + was 

30 stronger. 
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Table 7. Thiophene Adsorption Amount Normalized by Cation Content. 



(thiophene/ cation molar ratio) 



5 


Adsorbent 


Ag-Y 


Cu-Y 




Cations for 7i-complexation 


Ag + 


Cu + 


10 


Cation Content in Zeolite (wt%) 


38.2 


8.29 




Amount Adsorbed (molar ratio) 


0.42 


0.92 



15 The vapor-phase isotherms of benzene/thiophene on Cu-Y and Ag-Y 

were not completely reversible within a practical time scale (1 hours to 3 hours), 
because of relatively strong interaction. However, the bond energies are of the order of 
20 kcal/mol, which are weak enough so a substantial fraction of the adsorbed molecules 
are desorbed in a short period of time. For practical application, thermal desorption or 

20 displacement desorption may be a desirable option. 

Cu-X. In order to increase the Cu + content in the faujasite structure, an 
X-zeolite (Si/Al = 1.25) was ion-exchanged with Cu(N0 3 ) 2 . Cu-X was considered to be 
promising, since larger amounts of Cu + can be present in the structure of zeolite, 
especially at Sm sites. The color of Cu-X powder after ion-exchange (before auto- 

25 reduction) was bluer than that of Cu-Y, suggesting that more Cu 2+ was contained in X- 
zeolite than Cu-Y. For faujasite-type zeolites, the cations are designated as SI (the 
center of hexagonal prism), SI' (opposite SI but located in the cubooctahedron), SII 
(single six-ring in the supercage), SEP (opposite SII but inside the cubooctahedron), and 
SHI (near the four-ring windows in the supercage). See Yang, R. T. Gas Separation by 

30 Adsorption Processes, Butterworth, Boston (1987), Imperial College Press, London 
(1997); and Hutson, N. D.; Reisner, B. A.; Yang, R. T.; Toby, B. H., "Silver Ion- 
Exchanged Zeolites Y, X, and Low-Silica X: Observations of Thermally Induced 
Cation/Cluster Migration and the Resulting Effects on the Equilibrium Adsorption of 
Nitrogen," Chem. Mater. (2000), 12, 3020. The cations at SI, SI' and SIT sites are not 

35 exposed to the supercage and are shielded by framework oxygen. They cannot interact 
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with molecules inside the supercages directly. Only SII and SHI sites are exposed to the 
supercage and can interact with molecules inside the supercage. Because Y-zeolite 
(Si/Al=2.43) does not have the SIH cation sites, it is believed that only cations at SII 
site are responsible for the thiophene/benzene adsorption. On the other hand, X-zeolite 
5 (Si/Al=l .25) have a certain amount of SHI sites, based on the results on Ag-X. This 
additional Cu + located at SIH site should enhance the adsorption amount. However, 
unfortunately, it was found that Cu-X (Si/Al=1.25) was not stable at temperatures over 
200 °C, and the pore volume was decreased by 60 %. The color of the powder was 
changed from light blue to dark greenish gray after heating at 200 °C for 1 hour in 

10 vacuum. As a result, thiophene adsorption amount on Cu-X after 450°C auto-reduction 
was considerably smaller than that on Cu-Y (Isotherms were not shown). To alleviate 
the destruction of Cu-X structure, the heating rate to 450°C was reduced to 0.5 °C/min 
and 2°C/min (instead of 10°C/min ), but little improvement was observed. Although 
the reason for poor stability is not clear yet, copper cations were responsible for this 

1 5 phenomenon, because Ag-X, which will be discussed next, did not show any 

degradation of pore structure. It was reported that Nickel exchanged A-zeolite was 
unstable and lost its crystal structure at a temperature above 70°C. See Breck. D.W., 
Zeolite Molecular Sieves: Structure, Chemistry and Use; Wiley: New York (1974). 
The thermal instability of nickel exchanged A-zeolite was interpreted in terms of a 

20 ligand-field effect. After dehydration octahedral coordination for the nickel cation is 

not favored in the zeolite structure. A similar phenomenon might be occurred in copper 
exchanged X-zeolite. 

Ag-X. Thiophene adsorption performance on Ag-X was also examined 
to understand the effect of additional Ag + at SHI sites. Thiophene adsorption amounts 

25 on Ag-Y and Ag-X at 120°C were compared in Figure 8. Contrary to our expectation, 
Ag-X adsorbed less thiophene than Ag-Y on a per weight basis. However, since Ag-Y 
and Ag-X have different densities due to their different compositions, the adsorption 
amounts of thiophene were nearly the same based on per unit cell of faujasite (i.e., 
about 25 thiophene molecules/u.c. at 2.3x1 0" 5 atm and 36 - 40 thiophene molecules/u.c 

30 at 1 .8xl0" 2 atm for both Ag-Y and Ag-X). The reason for the ineffectiveness of extra 
Ag + at SHI site in Ag-X is not known. The thiophene adsorption amounts were also 
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examined from the viewpoint of zeolite pore volume and liquid thiophene volume. The 
molar density of liquid thiophene is calculated to be 83 cc/mol at 120°C. See Doong, 
S. J.; Yang, R. T., "A Simple Potential-Theory Model for Predicting Mixed-Gas 
Adsorption," Ind. Eng. Chem. Res. (1988), 27, 630. 1.5 mmol thiophene adsorbed per 
5 gram of Ag-Y at 2.3x1 0" 5 atm corresponded to 0.125 cc liquid thiophene per gram of 

Ag-Y. The thiophene liquid volume of 0.125 cc was nearly one half of pore volume of 
Ag-Y (0.24cc/g as shown in Figure 4), suggesting that one half of the pore was filled 
with thiophene. Therefore, it does seem likely that the cavity of Ag-X was filled with 
thiophene even at 2.3x1 0' 5 atm and only a small amount could be further adsorbed at 

10 higher pressures. One possible reason for the relatively low amounts in Ag-X is that 
thiophene molecules adsorbed could be located at positions near more than one silver 
cations at both SII and SHI sites simultaneously, so that no improvement in adsorbed 
amount was observed. 

Silver in Ag-X is known to form silver cluster at elevated temperature. 

15 See Sun, T. and K. Seff, "Silver Clusters and Chemistry in Zeolites," Chem. Rev. 

(1994), 94, 857. Hutson et al. showed that nitrogen adsorption capacity on Ag ion- 
exchanged low-silica X-zeolite (Si/Al = 1) was enhanced by heating in vacuum at 
450 °C for more than 4 hours. This enhancement was explained by thermally induced 
Ag cation and/or cluster migration using Rietvelt refinement of neutron powder 

20 diffraction data. The thiophene adsorption isotherm on Ag-X(Si/Al = 1 .25) after 450 °C 
for 4 hours is also shown in Figure 5. No enhancement was observed for thiophene 
adsorption on Ag-X(Si/Al = 1.25). 

H-USY. Thiophene and benzene isotherms on H-USY(Si/Al = 195) are 
shown in Figure 9. The interaction between thiophene/benzene and high-silica H-USY 

25 was very weak. Neither thiophene nor benzene was adsorbed at pressures lower than 
lxl 0" 3 atm. Adsorbed amounts were substantially smaller than Na-Y even at 2x10 2 
atm. These results are clearly due to the strong interactions between the cations in 
zeolite and the thiophene or benzene molecule. 

Na-ZSM-5. As stated above, several groups used Na-ZSM-5 for 

30 thiophene/benzene purification and showed that Na-ZSM-5 could remove thiophene 

impurities from benzene in their fixed bed breakthrough experiments. Figure 10 shows 
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isotherms on Na-ZSM-5. Although thiophene and benzene isotherms turned out to be 
virtually the same on Na-ZSM-5, it may be important to note that thiophene isotherms 
were fairly flat over the pressure range in the application of purification. This small 
difference of adsorbed amounts at 3xlO" 5 atm and 2x10 2 atm was preferable for the 
5 purification of benzene by removal of thiophene, which will be discussed below. One 
disadvantage of Na-ZSM-5 was small adsorption capacity owing to the small pore 
volume, as clearly seen in Figure 4. 

Activated carbon and modified activated alumina. The thiophene and 
benzene adsorption isotherms on activated carbon (Type PCB) and modified activated 

10 alumina (Selexsorb CDX) are shown in Figure 1 1 and Figure 12. PCB activated carbon 
is commercially designated for use in liquid phase or vapor phase applications such as 
recovery of alcohols, chlorinated hydrocarbons, hydrocarbons and aromatics. Selexsorb 
CDX is specially formulated by Alcoa Industrial Chemicals for adsorption of polar 
organic compounds including sulfur-based molecules (mercaptans, sulfides, disulfides, 

15 thiophenes), nitrogen-based molecules (nitriles, amines, pyridines) and oxygenated 

hydrocarbon molecules (alcohol, glycols, aldehydes, ketons, ethers, peroxides). More 
benzene than thiophene was adsorbed on both activated carbon and modified activated 
alumina at pressures below 1x10 2 atm. This is surprising because Selexsorb was 
thought to adsorb thiophene more selectively than benzene. This result indicates that 

20 polarizability may play a more important role (than permanent dipole moment and 

magnetic susceptibility) in adsorption on these two sorbents. At pressures higher than 1 
x 10" 2 atm, pore filling dominated hence benzene adsorbed less than thiophene because 
of the larger molar volume of benzene than thiophene. 

Comparison of thiophene adsorption on all sorbents. Thiophene 

25 adsorption isotherms on all sorbents are compared in Figure 13. It is clearly seen that 
Ag-Y and Cu-Y could adsorb significantly larger amounts of thiophene even at very 
low pressures. 

Purification capability 

30 For purification, as has been shown for the removal of ppm levels of 

dienes from olefins, an important consideration is the separation factor; and the 
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separation factor should be evaluated at low concentrations of the solute but at high 
concentrations of the solvent. In other words, the capacity for the dilute solute (at ppm 
level) is a factor that determines the purification. As stated above, the separation factor 
(a) is given by: 

5 

a n = yi x Y 2 /{X 2 Y x ) (3) 

where X and Y are mole fractions in the adsorbed and fluid phases, respectively. In our 
case, component 1 is thiophene and component 2 is benzene. A further important 
10 parameter is the capacity of the sorbent for the solute at low concentrations. Hence, the 
ability of the sorbent for purification depends on the amount of thiophene adsorbed at 
the ppm levels. 

To estimate the adsorption amount in the mixture of thiophene and 
benzene, the pure-component adsorption data were first fitted by both D-A (Dubinin- 

1 5 Astakhov) and L-F (Langmuir-Freundlich) isotherms. The separation factors of 
thiophene over benzene were then obtained using the equivalent multi-component 
isotherms, as described previously. The conversion of molar adsorbed amount to 
volumetric adsorbed amount and vice versa were performed in the same way as Doong 
and Yang. The molar volumes of thiophene in this work were 79.7 cc/mol at 90°C, 

20 83.0 cc/mol at 120°C and 89.8 cc/mol at 180°C, and those of benzene were 90.5 

cc/mol at 90°C, 94.6 cc/mol at 120°C and 102.9 cc/mol at 180°C. The adsorbed molar 
volumes were assumed to be the same in both single- and mixed-gas adsorption. 
Separation factors were calculated for two cases: lxl0" 4 atm thiophene/1 atm benzene 
and lxl 0" 5 atm thiophene/1 atm benzene. 

25 The parameters for the D-A and L-F equations are given in Table 8 and 

the resulting separation factors are listed in Table 9. Some of the D-A parameters were 
not applicable for the Doong/Y ang model for mixtures, when the volumetric adsorbed 
amount of benzene was much higher than that of thiophene. The separation factors of 
thiophene over benzene follow the order: Ag-Y > Na-ZSM-5 > Cu-Y « activated 

30 carbon » Na-Y » H-USY « modified activated alumina. This order can be explained 
as follows. As described earlier, the ability of the sorbent for purification depends on 
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the amount of thiophene adsorbed at the ppm levels, compared with the amount of 
benzene adsorbed at high concentration. Ag-Y and Na-ZSM-5 showed fairly flat 
isotherms of thiophene over the entire partial pressure range, resulting in the higher 
separation factors. In fact, Ag-Y and Na-ZSM-5 maintained 65% and 47% of 
thiophene adsorption amounts at 2.3xl0~ 5 atm compared with those at 1.3xl0" 2 atm, 
respectively. The advantage of Ag-Y against Na-ZSM-5 was further enhanced by the 
fact that Ag-Y was slightly selective for thiophene than benzene, while 
thiophene/benzene isotherms on Na-ZSM-5 were virtually the same. On the other 
hand, Cu-Y and activated carbon could adsorb a certain amount of thiophene/benzene, 
but smaller adsorption amounts at lower pressures reduced the separation factors 
considerably. Only 18% and 5% of thiophene adsorption amounts were retained by Cu- 
Y and activated carbon between 2.5 x 10 5 atm and 1.3x10 2 atm, respectively. Larger 
adsorption amounts of benzene than thiophene by Selexsorb CDX made the separation 
factor of this sorbent extremely low. 
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Table 8. Parameters for the D-A and L-F Equations 





L-F Equation 
qo B n 
(mmol/g) (1/atm) (-) 


1J-A equation 
qo Vo C n Ps 
(mmol/g) (cc/g) (-) (-) (atm) 


Ag-Y Benzene 120°C 
180°C 

Thiophene 120°C 
180°C 


3.133 3.40 7.67 
2.032 155.5 1.85 
3.200 4.416 6.59 
2.602 6.13 4.77 


2.416 0.2286 0.0559 1.56 2.96 
2.008 0.2066 0.0786 4.57 10.07 
2.986 0.248 0.0609 1.12 2.68 
2.329 0.2091 0.0699 1.87 9.42 


Cu-Y Benzene 90 °C 
120°C 
Thiophene 90 °C 
120°C 


2.989 292.6 1.76 
2.897 178.6 1.59 
3.910 196.4 1.80 
3.548 121.5 1.65 


2.940 0.266 0.0889 4.36 1.34 
2.929 0.2771 0.1008 3.95 2.96 
3.826 0.3048 0.0935 3.99 1.19 
3.548 0.2946 0.1010 3.81 2.68 


Na-Y Benzene 120°C 
180°C 
Thiophene 120°C 
180°C 


2.781 13985 0.89 
2.820 449.0 0.90 
2.681 13548 0.84 
2.739 15001 0.50 


2.789 0.2639 0.0980 6.80 2.96 
2.872 0.2956 0.1200 4.73 10.07 
2.658 0.2207 0.1050 7.30 2.68 
2.846 0.2555 0.1330 7.54 9.42 


PCB-AC Benzene 90°C 
120°C 

Thiophene 90 °C 
120°C 


7.828 4.47 2.46 
6.312 4.254 2.16 
7.231 5.54 2.29 
6.996 6.39 1.66 


6.329 0.5725 0.1638 1.71 1.34 
6.197 0.5863 0.1847 1.62 2.96 
8.936 0.7118 0.2246 1.24 1.19 
7.236 0.6009 0.2193 1.59 2.68 


Na-ZSM5 Benzene 90 °C 
120°C 

Thiophene 90 °C 
120°C 


0.928 17.15 3.43 
0.795 14.66 3.47 
1.023 16.91 3.27 
0.811 16.36 3.30 


0.987 0.0893 0.0810 1.44 1.34 
0.826 0.0781 0.0793 1.64 2.96 
1.036 0.0825 0.0880 1.79 1.19 
0.856 0.0711 0.0832 1.70 2.68 


SelexsorbCDX Benzene 90 °C 
120°C 

Thiophene 90°C 
120°C 


1.182 1196.6 1.39 
1.096 75902 0.71 
1.223 29720 0.79 
1.134 133016 0.55 


1.176 0.1064 0.0914 5.23 1.34 
1.109 0.1049 0.1049 8.14 2.96 
1.239 0.0987 0.1141 6.78 1.19 
1.158 0.0962 0.1265 8.20 2.68 


H-USY Benzene 90°C 
120°C 

Thiophene 90 °C 
120°C 


4.258 568.1 0.48 
3.895 63.573 0.53 
4.568 382.1 0.51 
4.356 24.41 0.61 


4.359 0.3943 0.2800 4.21 1.34 
4.100 0.3879 0.2825 3.63 2.96 
4.689 0.3735 0.2939 3.55 1.19 
4.599 0.3819 0.2820 4.07 2.68 
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Table 9. Separation Factors estimated by the D-A and L-F Equations at 120° C. 



Thiophene(atm)/B enzene(atm) 


L-F Equation 
1E-4/1 1E-5/1 


D-A Equation 
1E-4/1 1E-5/1 


Ag-Y 


3300 


23000 


1600 


11000 


Cu-Y 


31 


78 


360 


700 


Na-Y 


0.16 


0.10 


*** i) 


*** i) 


PCB-AC 


64 


160 


33 


83 


Na-ZSM5 


700 


3500 


*** i) 


*** i) 


Selexsorb CDX 


9.6E-04 


1.5E-04 


*** i) 


*** 0 


H-USY 


1.2E-03 


2.7E-04 


*** i) 


*** i) 



Not application to Doong/Y ang extension. 



Another important parameter for purification is the sorbent capacity for the 
solute at very low concentrations, i.e., for thiophene at ppm levels. The sorbent 
capacities are shown in Figure 13. At 2.3 xlO 5 atm, the capacities for thiophene follow 
the order: Ag-Y & Cu-Y » Na-ZSM-5 > activated carbon > Na-Y > modified alumina 
& H-USY. Based on these two key parameters, the 7t-complexation sorbents are clearly 
the best. 

Heat of adsorption 

Heats of adsorption were calculated using the Clausius-Clapeyron 
equation from isotherms at two different temperatures, and are shown in Table 10. All 
the heats of thiophene/benzene adsorption had the tendency to decrease as the loading 
increased. This is a common phenomenon for the sorbents such as ion-exchanged 
zeolites that have heterogeneous sites. The heats of adsorption on activated carbon, in 
particular, ranged widely from 23.9 kcal/mol (at 0.5mmol/g loading) to 8.0 kcal/mol (at 
3mmol/g). Ag-Y and Cu-Y exhibited the higher heats of adsorption than Na-Y for both 
benzene and thiophene, it is believed due to rc-complexation. More importantly, the 
heats of adsorption for thiophene were higher than that of benzene. These experimental 
results can be explained by molecular orbital calculation and NBO analysis, which will 
be discussed shortly. 

At the low loading of 0.5 mmol/g, Na-ZSM-5 showed nearly the same 
heats of adsorption as Na-Y for both thiophene and benzene. The different pore 
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dimensions (5.2-5.6 A for ZSM-5 vs. 7.4 A for Na-Y) apparently had no influence on 
the heats of adsorption. It is not clear why the amounts adsorbed on Na-Y decreased 
sharply at very low pressures while that on Na-ZSM-5 maintained. 



Table 10. Heat of adsorption (kcal/mol) calculated from isotherms at different temperatures. 





Na-Y 


Ag-Y 


Cu-Y 


H-USY 


Na-ZSM5 


Activated Carbon Activated Alumina 




(Si/Al=2.43) (Si/Al=2.43) (Si/Al=2.43) (Si/Al=195) (Si/Al=10) 


Type PCB 


Selexsorb CDX 


Thiophene 


19.1-19.6 


21.3-21.5 


20.8-22.4 


7.9-11.2 


18.6-19.2 


8.0-23.9 


16.1-17.5 




(0.5-2.0) 


(1.5-1.7) 


(2.0-3.0) 


(0.1-0.3) 


(0.45-0.60) 


(0.5-3.0) 


(0.2-1.0) 


Benzene 


17.0-18.2 


19.0-20.1 


19.3-21.8 


6.6-13.1 


16.5-17.9 


13.1-16.1 


16.8-19.6 




(1.5-2.0) 


(1.5-1.8) 


(1.8-2.5) 


(0.1-0.5) 


(0.45-0.65) 


(1.0-3.0) 


(0.6-1.0) 



*) Numbers in parentheses indicate the adsorption amounts (mmol/g) for calculation. 
**) The heats of adsorption decreased with loading in all cases. 



Bond Energies, Geometries and NBO results 

The energies of adsorption calculated using equation 1 for thiophene and 
benzene are summarized in Table 1 1 . The theoretical calculations indicate that the n- 
complexation strengths follow the order CuZ > AgZ and more importantly, thiophene > 
benzene. This trend is in agreement with the experimental data, in Table 10. In fact, 
the molecular orbital results on CuZ and AgZ are in excellent agreement with the 
experimental data. Both chloride and zeolite models were used as the anion in the 
theoretical calculations, while only zeolite framework was the anion in the experiment. 
It is believed that the anion has a large effect on the 7t-complexation bonds. The bond 
energies on the zeolites (Z) are significantly higher than that on the chlorides (Table 
11). This result indicates that the zeolite anion is more electronegative than the chloride 
anion, which has been revealed by Chen and Yang in their ab initio molecular orbital 
calculations. In the optimized structures of thiophene-MCl complexes, the distance 
between the thiophene molecule and Cu ion is about 0.3 A shorter than that of 
thiophene and Ag ion for chloride. The NBO analysis is summarized in Tables 12 and 
13. There is some donation of electron charges from the 7i-orbital of thiophene to the 
vacant s orbital of metals known as a donation and, simultaneously, back-donation of 
electron charges from the d orbitals of metals to n orbital (i.e., anti-bonding n orbital) 
of thiophene or n back-donation. It appears that the a donation is more predominant for 
thiophene and the it back-donation is more important for benzene. Comparing the two 
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anions, zeolite anion and chloride anion, the NBO results show that both a donation and 
d - n back-donation are significantly stronger with the zeolite anion bonded to Ag + or 
Cu + . The charge transfer results again confirmed the experimental data that the relative 
strengths of the 7C-complexation bonds follow the order: thiophene > benzene and Cu + > 
Ag + . 

Table 11. Summary of energies of adsorption for thiophene and benzene in kcal/mol calculated 



from molecular orbital theory (Z denotes Zeolite anion) 




E^CThiophene) 


E^Benzene) 


CuCl 


13.5 


12.4 


AgCl 


9.0 


8.6 


CuZ 


21.4 


20.5 


AgZ 


20.0 


19.1 


Table 12. Summary of NBO analysis* of 7r-compIexation between thiophene and MC1/MZ 




C-M interaction 


M-C interaction 




(a donation) 


(d - n* back-donation) Net Change 


MX 


Qi 


<L qi+<l2 


CuCl 


0.037 


-0.022 0.015 


AgCl 


0.022 


-0.014 0.008 


CuZ 


0.112 


-0.063 0.049 


AgZ 


0.101 


-0.086 0.015 



* qj is the amount of electron population change in valence s orbitals of the metal, and is the total 
amount of electron population decrease in valence d orbitals of the metal. 



Table 13. Summary of NBO analysis* of rc-complexation between benzene and MC1/MZ. 





C-M interaction 


M-C interaction 






(a donation) 


(d - 7i* back-donation) 


Net Change 


MX 


qi 


<h 


Qi +02 


CuCl 


0.011 


-0.013 


-0.002 


AgCl 


0.003 


-0.007 


-0.004 


CuX 


0.078 


-0.119 


-0.041 


AgX 


0.043 


-0.045 


-0.002 



* q, is the amount of electron population change in valence s orbitals of the metal, and is the 
total amount of electron population decrease in valence d orbitals of the metal. 
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Uptake Rates 

Table 14 compares diffusion time constants for Ag-Y, Cu-Y and Na-Y. 
These constants were calculated from the uptake rates using the solution for Fick's 
law. Uptake rates among the three sorbents were very similar and no notable 
differences between thiophene and benzene were observed. These rates are reasonably 
fast for practical applications. 

In this experiment, vapor-phase benzene/thiophene adsorption 
isotherms were investigated to develop new sorbents for desulfurization. Among the 
sorbents studied, it was found that Cu-Y and Ag-Y exhibited excellent adsorption 
performance (capacities and separation factors) for desulfurization. This enhanced 
performance compared to Na-Y is believed to be due to the 7i-complexation with Cu + 
and Ag + . Molecular orbital calculations confirmed the relative strengths of 7t- 
complexation: thiophene > benzene and Cu + > Ag + . In an actual desulfurization 
process, liquid phase adsorptions such as batch adsorption or fixed-bed adsorption 
column around room temperatures are viable options. Liquid phase experiments are 
described herein. Even in the case of liquid phase separation, enhanced selectivities 
based on 7r-complexation studied in this experiment are applicable. 



Table 14. Diffusion Time Constants (1/s) of Thiophene and Benzene. 





Thiophene 




Benzene 






Pressure Change 


D/r 2 


Pressure Change 


D/r 2 




(atm) 


(1/s) 


(atm) 


(1/s) 


Ag-Y 


1.9E-5- 1.6E-4 


2.0E-05 


5.7E-5 - 5.5E-4 


5.4E-05 




5.6E-4- 1.8E-3 


3.3E-04 


5.5E-4 - 5.4E-3 


L8E-04 




1.8E-3- 1.5E-2 


1.2E-04 


5.4E-3 - 2.6E-2 


4.2E-04 1 


Cu-Y 


2.3E-5 - 8.8E-5 


1.3E-05 


3.2E-5- 1.2E-4 


2.0E-05 




8.8E-5 - 4.0E-4 


2.3E-05 


1.2E-4-5.4E-4 


5.2E-05 




4.0E-4 - 2.2E-3 


1.5E-04 


5.4E-4-3.0E-3 


3.3E-04 




2.2E-3 - 1.3E-2 


1.4E-04 


3.0E-3 - 1.8E-2 


3.7E-04 


Na-Y 


8.8E-5 - 4.0E-4 


1.7E-05 


0- 1.5E-4 


1.5E-05 




4.0E-4 - 2.2E-3 


1.6E-04 


1.5E-4-5.5E-4 


2.6E-05 




2.2E-3 - 1.3E-2 


2.8E-04 


5.5E-4 - 5.4E-3 


5.6E-05 




1.3E-2-2.9E-2 


5.4E-04 


5.4E-3 - 2.6E-2 


2.4E-04 



EXPERIMENT D 

As mentioned hereinabove, it has further been fortuitously found that the 
sorbents of the present invention also are excellent sorbents for removal of aromatics. A 
5 detailed description of removal of aromatics follows. 

Adsorption of benzene and cyclohexane on various Y-zeolites was 
investigated in order to develop sorbents for the purification of aliphatics by removal of 
aromatics. Ag-Y showed superior benzene/cyclohexane selectivities than Na-Y, Pd-Y 
and H-USY. Separation factors greater than 10 5 were obtained with Ag-Y at low 

10 concentrations of benzene. The high selectivities were achieved by the strong 

interaction between benzene and Ag-Y, while the interaction with cyclohexane was not 
influenced by the cations. Molecular orbital calculation revealed that benzene formed 
classical 7i-complexation bond with Ag-Y: donation of electron charges from the p- 
orbitals of benzene to the vacant s-orbital of the silver (a donation) and, simultaneously, 

15 back- donation of electron charges from d-orbitals of silver to 7t* orbital of benzene (d- 
7i* back- donation). Grand canonical Monte Carlo simulations were also performed for 
the adsorption isotherms. Potential parameters of benzene on Ag-Y including iz- 
complexation were first developed. GCMC simulation and experiments were in 
excellent agreement for adsorption of benzene on Y-zeolites. 

20 Purification of aliphatics by removing aromatics is important in both 

petrochemical industry and for pollution control. Because benzene (and aromatics) is 
known to be highly carcinogenic, concentration of benzene needs to be minimized in the 
automotive fuel. A number of separation processes have been employed to reduce the 
benzene concentration in a refinery's gasoline pool, so that it meets new reformulated 

25 gasoline requirements. In a typical benzene reduction process, a combination of 

extraction and distillation is used (Jeanneret, 1997). Also, removal of aromatics from 
kerosene improves the clean burning properties of the fuel and separation of the 
aromatics from the isoparaffins and naphthenes in lubricating oils improves the 
viscosity-temperature relationship (Bailes, 1983). 

30 Because of the importance of aromatics/aliphatics separation and the 

problems associated with current separation processes, possible alternatives have been 
under continuing investigation. These include pervaporation (Hao et al, 1997), liquid 
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membranes (Li, 1968), and the use of liquid inclusion complexes (Atwood, 1984). 
Purification of dilute aromatics from aliphatics (e.g., toluene and/or xylene in heptane) 
by temperature swing adsorption (TSA) was studied in the liquid phase (Matz and 
Rnaebel, 1990). Commercially available sorbents were used: silica gel, activated 
5 alumina, activated carbon, zeolite 13X, and polymeric resin (XAD-7). Among these 
sorbents, silica gel was considered the best due to its superior thermal-exchange 
capacity. However, the selectivities were low. 

Adsorption is playing an increasingly important role in separations. 
However, its utility is limited by the availability of selective sorbents. The conventional 

10 sorbents and separations using them are based on van der Waals and electrostatic 

interactions between the sorbate and the sorbent. Studies of new sorbents using weak 
chemical bonds such as chemical complexation have begun only recently in our 
laboratory. As suggested by King (King, 1987), chemical complexation bonds are 
generally stronger than van der Waals interactions, yet weak enough to be reversible. 

15 Therefore, tremendous opportunities exist for developing new sorbents and new 

applications in separations by using weak chemical bonds, including various forms of 
complexation bonds. 

The 7r-complexation is a sub-class of chemical complexation. It pertains 
to the main group (or d-block) transition metals, i.e., from Sc to Cu, Y to Ag, and La to 

20 Au in the periodic table (Cotton and Wilkinson, 1966). These metals or their ions can 
form the normal cr-bond to carbon and, in addition, the unique characteristics of the d 
orbitals in these metals or ions can form bonds with unsaturated hydrocarbons in a 
nonclassical manner. This type of bonding is broadly referred to as 7c-complexation. 
This 7r-complexation has been seriously considered for olefin/paraffin separations using 

25 liquid solutions containing silver or cuprous ions (Quinn, 1971; Ho et al., 1988; Keller 
et al., 1992; Blytas, 1992; Eldridge, 1993; Safarik and Eldridge, 1988). These involved 
gas-liquid operations. While gas-solid operation can be simpler as well as more 
efficient, particularly by pressure swing adsorption, the list of attempts for developing 
solid Tc-complexation sorbent is a short one. 

30 More recently, several new sorbents based on 7i-complexation were 

prepared for selective olefin adsorption. These included Ag + -exchanged resins (Yang 
and Kikkinides, 1995; Wu et al., 1997), monolayer CuCl on pillared clays (Cheng and 
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Yang, 1995) and monolayer AgN0 3 /Si0 2 (Padin and Yang, 1997; Rege et al., 1998; 
Padin and Yang, 2000). Moreover, Ag ion-exchanged Y-Zeolite was developed to 
purify butene by removing trace amounts of butadiene (Padin et al., 1999). The 
purification performance of this sorbent was superior to that of NaY owing to n- 
5 complexation. 

It is possible, in principle, to purify aromatics from aliphatics based on 
7C-complexation. In the benzene molecules, the carbon atom is sp 2 hybridized. Hence, 
each carbon has three sp 2 orbitals and another P z orbital. The six P z orbitals in the 
benzene ring form the conjugative n bonds. The P z orbitals also form the antibonding 

10 7i* orbitals, which are not occupied. When benzene interacts with transition metals, the 
7r-orbitals of benzene could overlap with the empty or unfilled outer-shell s orbital of the 
transition metal to form a a-bond. Furthermore, the vacant antibonding 7t*-orbitals of 
benzene could overlap with the d-orbitals in the transition metals, similar to that formed 
in the olefin-Cu + bond (Huang et al., 1999a). In this work, Ag + and Pd 2+ ion-exchanged 

15 zeolites were selected for the study of aromatics/aliphatics separations, since these two 
cations have the most promising 7i-complexation capabilities, and they are stable. In 
order to understand the effects of Tc-complexation, Ag + and Pd 2+ ion-exchanged zeolites 
were compared with that of Na + ion-exchanged zeolite and a high-Si0 2 zeolite. . 

Benzene and cyclohexane are an ideal pair of model compounds for 

20 studying selective sorbents for purification. These molecules have similar shapes, 

polarizabilities and volatilities (the boiling points are 80°C for benzene and 81 °C for 
cyclohexane). The kinetic diameter of benzene, which is calculated from the minimum 
equilibrium cross-sectional diameter, is estimated to be 5.85 angstrom compared with 
6.0 angstrom for cyclohexane (Breck, 1974). Therefore, benzene and cyclohexane were 

25 used in this work. 

In addition to the experimental investigation, molecular orbital 
calculations were performed to obtain a basic understanding of the origin of the strong 
interactions of Ag ion-exchanged zeolite with benzene. Moreover, grand canonical 
Monte Carlo (GCMC) simulations were used to simulate/predict the isotherms. Many 

30 GCMC studies have been performed for predicting adsorption isotherms in zeolites. 

These include N 2 , 0 2 , Ar, Xe, CH 4 , C 2 H 4 and C 4 H 10 on zeolite- A, X, Y and MFI-type 
Zeolite (Woods and Rowlinson, 1989; Razmus and Hall, 1991; Kravias and Myers, 
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1991; Watanabe et al., 1995). Adsorption of benzene on NaY, Heulandite and MFI-type 
zeolites has also been studied (Snurr et al, 1994; Kitagawa et al., 1996; Klemm et al., 
1998; Laboy et al., 1999). However, these studies did not involve the effects of k- 
complexation on adsorption. In this study, benzene adsorption isotherms on Ag ion- 
5 exchanged Y-zeolite were simulated. 

Experimental Details 

Sorbent Preparation 

The sorbents in this work were ion-exchanged zeolites. From earlier 

10 results on olefin/paraffin separations (Yang and Kikkinides, 1995; Wu et al., 1997; 
Cheng and Yang, 1995; Padin and Yang, 1997; Rege et al., 1998; Padin and Yang, 
2000), Ag + and Pd 2+ were the most promising ions because of their strong 71- 
complexation abilities. 

Ion-exchange of zeolites was accomplished with aqueous solutions. 

15 Faujasite was selected as the zeolite framework structure for ion-exchange, because of 
its large pore aperture (7.4 angstrom), large pore volume (0.489 ml/ml) and wide range 
of Si0 2 /Al 2 0 3 ratios(l-oo). Two as-received zeolites, Na-type Y-Zeolite (Na-Y, 
Si/Al=2.43, Strem Chemical) and H-type ultra-stable Y-Zeolite (H-USY, Si/Al=195, 
TOSOH CORPORATION, HSZ-390HUA), were used in this study. Pd 2+ type Y-Zeolite 

20 was prepared by exchanging Na-Y zeolite with Pd(N0 3 ) 2 (Aldrich) in an aqueous 

solution (Li et al., 1994). First, Na-Y was exchanged using 2-fold excess of Pd(N0 3 ) 2 
(0.005M) at 80 °C for 24 hours. After the exchange, the zeolite suspension was filtered 
and washed with copious amount of deionized water. The product was dried at 1 10°C 
overnight in air. Ag + ion-exchange was prepared at room temperature for 24 hours in 

25 the same manner as Pd 2+ exchange, using 2-fold excess AgN0 3 (0.1M). 

Adsorption Isotherms and Uptake Rates 

Single component isotherms for benzene and cyclohexane were measured 
using standard gravimetric methods. A Shimadzu TGA-50 automatic recording 
30 microbalance was employed. Helium (Pre-purified grade, Metro Welding, 99.995%) 
was used as the carrier gas and was first passed through two consecutive gas-wash 
bottles (to ensure saturation) which contained benzene (HPLC grade, Sigma- Aldrich, 
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99.9+%) or cyclohexane (HPLC grade, Aldrich, 99.9+%). After diluting the 
concentration to the desired value by blending with additional helium (total flowrate: 
250 cc/min), the mixture was directed into the microbalance. 

Isosteric heats of adsorption were calculated using the Clausius- 
5 Clapeyron equation from isotherms at different temperatures. Nitrogen isotherms at 
77K measured with a Micromeritics ASAP 2010 system were used for pore size 
distribution and pore volume determination. Pore size distributions were calculated 
with Horvath-Kawazoe equation (Horvath and Kawazoe, 1983; Cheng and Yang, 1994). 
Pore volumes were calculated at 0.95 P/Po. The overall diffusion time constants, D/r 2 , 
10 were calculated from the uptake rates (Yeh, 1989). In this work, the short time region 
was used. 



Molecular Orbital Calculation Details 

The 7t-complexation bonding between olefins and sorbents have been 
1 5 extensively investigated using molecular orbital (MO) calculations (Chen and Yang, 
1996; Huang et al., 1999a, 1999b). More recently, MO studies were extended to 
benzene adsorption on various salts such as CuCl and AgCl (Takahashi et al.). In this 
work, similar MO studies for benzene adsorption on Ag ion-exchanged zeolite were 
performed using a cluster model of Ag-zeolite. The Gaussian 94 Progam (Frisch et al., 
20 1995) in Cerius2 molecular modeling software (Bowie, 1995) from Molecular 

Simulation, Inc. was used for calculations. MO calculations for benzene and Ag zeolite 
were performed at the Hartree-Fock (HF) and density functional theory (DFT) level. 
The 3-21G basis set is the split- valence basis set, which has been used successfully for 
many simulations for benzene adsorption on zeolites (O'Malley and Braithwaite, 1995), 
25 protonation reaction of propylene and isobutene on zeolite (Viruera-Martin et al., 1993) 
and 0 2 , N 2 and C 2 H 4 adsorption on Ag-Zeolites (Chen and Yang, 1996). The reliability 
of this basis set has been confirmed by the accuracy of calculation results in comparison 
with experimental data. Therefore, the 3-2 1G basis set was employed for both geometry 
optimization and natural bond orbital (NBO) analysis. 



30 
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Geometry Optimization and Energy of Adsorption Calculations 

The restricted Hartree-Fock (RHF) theory at the 3-21G level basis set 
was used to determine the geometries and the bonding energies of benzene on Ag zeolite 
cluster, since Ag + has filled d-orbitals with spin = 1 . 
5 The computational model used for Ag zeolite were the same as that used 

by Chen and Yang (Chen and Yang, 1996), since this model yielded heats of adsorption 
of N 2 and C 2 H 4 on Ag zeolite that agreed well with experimental data. The simplest 
cluster model, (HO) 3 Si-0(Ag)-Al(OH) 3 , was chosen for this work. The optimized 
structures were then used for bond energy calculations according to the following 
10 expression: 

F = F 4- F - F (1) 

^ads -^adsorbate ^adsorbent ~ -^adsorbent-adsorbat 

Where E adsorbate is the total energy of benzene, E adsorbent is the total energy of the bare 
15 adsorbent i.e., the Ag zeolite cluster, and 

Eadsorbem-adsorbate is the total energy of the 

adsorbate/adsorbent system. 

Natural Bond Orbital (NBO) Analysis 

The optimized structures were also used for NBO analysis at the 

20 B3LYP/3-21G level. The B3LYP (Becke, 1993a) approach is one of the most useful 

self-consistent hybrid (SCH) approaches (Becke, 1993b), which is Becke's 3-parameter 
nonlocal exchange functional (Becke, 1992) with the nonlocal correlation functional of 
Lee, Yang and Parr (Lee et al., 1988). 

The NBO analysis performs population analysis pertaining to the 

25 localized wave-function properties. It gives a better description of the electron 

distribution in compounds of high ionic character, such as those containing metal atoms 
(Reed et al., 1985). It is known to be sensitive for calculating localized weak 
interactions, such as charge transfer, hydrogen bonding and weak chemisorption. 
Therefore, the NBO program (Glendening et al., 1995) was used for studying the 

30 electron density distribution of the adsorption system. 



49 

Monte Carlo Simulation Details 

Grand canonical ensemble (at a fixed pressure) Monte Carlo simulation 
was performed utilizing the Cerius2 molecular modeling software (Bowie, 1995) on a 
Silicon Graphics Indigo2 workstation running IRIX v.6.5. After the "crystal builder" 
5 module was used to set up the zeolite models, energy expressions and parameters for 
force field were input by the "force field editor" module. The "sorption" module was 
then used for benzene adsorption simulations. 

The simulation was performed at 120°C and 180°C for at least 1,000,000 
configurations. During the simulation, the sorbate molecules interacted with the 

10 potential field generated by sorbate-sorbent and sorbate-sorbate interactions. The initial 
configurations contained no sorbate molecules. Each subsequent configuration was 
generated by one of four moves: CREATE, DESTROY, TRANSLATE and ROTATE. 
Molecular creation attempts were made at random points within the accessible portion 
of the zeolite lattice with the criterion that there be no overlapping sites (i.e., creations 

1 5 which results in interaction sites which were closer than half the sum of the van der 
Waal's radii of the two sites were rejected). Translation and rotation were performed 
with the same rejection criterion. For non-overlapping sites, the change in the potential 
energy accompanying the new configuration was calculated and subsequently accepted 
or rejected in accordance with standard acceptance probabilities (Molecular Simulation 

20 Inc., 1998). 

Zeolite Models 

The structure determined from powder neutron diffraction data by Fitch 
et al. (Fitch et al., 1986) was used as the structure for Na-Y Zeolite (Si/Al=2.43). The 

25 location of the Na cations were 32 atoms on SII sites, 16 atoms on SI sites and 8 atoms 
on SF sites. Na cations on SI and SF sites were located in such a manner that two 
adjacent SI and SF were not simultaneously occupied. This is due to the electrostatic 
repulsion caused by their small separation of 0.218 nm (Fitch et al., 1986). The position 
of Al in the framework was randomly chosen using the "disorder" function in the 

30 "crystal builder" module that obeyed the Loewenstein rule. No charge-free "blocking 
atom" was added because large molecules such as benzene were not located in the site, 
which was actually sterically inaccessible. 
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For the simulation of Ag ion-exchanged zeolite (Si/Al=2.43), the Na 
cations in Na-Y were replaced by Ag maintaining the same cation locations. An all silica 
Y-Zeolite structure was used to simulate Ultra-Stable Y-Zeolite (Si/Al=195) since this 
zeolite has less than one Al atom per unit cell. 

Until recently, the location of the Ag cation was not accurately 
determined. Using neutron diffraction data, Hutson et al (Hutson et al., 2000) reported 
the Ag location and occupancies of Ag in Ag-Y (Si/Al=2.43) after 450°C dehydration. 
Their results showed that the Ag cation sites were 28 atoms on SII sites, 4 atoms on SIT 
site, 1 1 atoms on SI sites and 12 atoms on SF sites. Simultaneous occupancy of Ag 
cations at SII and SII' as well as SI and SF sites was unlikely due to the large repulsion. 
This location is very similar to the location of Na cations in Na-Y reported by Fitch et 
al., 1986. Therefore, the zeolite model for Ag-Y used in this work is considered to 
reflect the actual Ag cation locations in Ag-Y (Si/Al=2.43). 

Forcefield Determination 

One of the most important factors for GCMC simulation is the selection 
of an energy expression and parameters for the forcefield. The forcefield used in these 
simulations was a modified version of Cerius2 Sorption Demontis Forcefield, since 
Demontis et al. simulated the mobility of benzene in Na-Y zeolite successfully using 
molecular dynamics methods (Demontis et al., 1989). In this model, the total potential 
energy (U Total ) between the zeolite lattice and adsorbate molecules is written as the sum 
of the interactions between adsorbate molecules(AA) and that between the adsorbate 
molecules and the zeolite lattice (AZ). 

Urotai- Uaa+ Uaz (2) 

Interactions between adsorbate molecules (Uaa) were described as the sum of the atom- 
atom Buckingham potentials (the first and second terms below), which have shown to 
yield reasonable predictions for liquid and solid benzene (Williams and Cox, 1984), and 
the Coulombic interaction potentials (the third term below). 
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C„ 



(3) 



The first and second terms represent, respectively, the repulsive and dispersive potential 
energies between sites "i" and "j", and the third term represents the interaction potential 
5 between point charges q { and qj of sites "i" and "j" separated by a distance r^. 

The interactions between the adsorbate molecules and the zeolite lattice 
(Uaz) * s written as the sum of the atom-atom Lennard- Jones 12-6 potential and the 
Coulombic interaction potentials. 
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10 

Following the empirical approach used by Razman and Hall (Razmus and 
Hall, 1991), Watanabe et al. (Watanabe et al., 1995) used the empirical approach to 
parameterize the forcefield for adsorption of N 2 , 0 2 , and Ar in type A and faujasite 
zeolites. The dispersive part of the Lennard- Jones 12-6 potential (U dispersion ) was written 
1 5 using an adjustable parameter (P { ) of the adsorbate and the polarizabilities (a) of the 
atoms of the adsorbate molecules (i) and in the zeolite lattice (j). 



U 



dispersion 



r \ 
a i a j 

\ v J 



(5) 



The potential energy well depth parameter, is then written as 



*„ = Pi 



Act 



(6) 



v J 



20 The Lennard- Jones parameter, ay, for the adsorbate-zeolite interaction is given by the 
mixing rule, 



o-, = 



(7) 
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The Oj for the zeolite lattice are related to the van der Waals radii, Rj, of the respective 
ions by 

R j = 2^a j (8) 

5 

The well depth parameter, Cy , is assumed to be a geometric combination of e { and e j5 

) 05 (9) 

10 The van der Waals energy in the periodic framework was calculated using a minimum 

image convention with an interaction cut-off distance of 1 nm. The Coulombic term was 
evaluated using the Evald summation method (Allen and Tildesley, 1987). 

The potential parameters used for Na-Y, Ag-Y and H-US Y are listed in 
Table 15. The potential parameters for Na-Y are the same as that of Demontis et al 

15 (Demontis et al., 1989). Van der Waals interactions between adsorbate and zeolite cage 
were modeled by means of interactions between the (C,H) and (Na,0) atoms. 
Interactions with the Si and Al atoms were neglected because they were well shielded by 
the oxygen atoms of the Si0 4 and A10 4 tetrahedra. Electrostatic interactions between all 
charges were included in the potential energy calculation. As for Ag-Y and all silica Y 

20 zeolites, modification of parameters was performed using the mixing rule and the 

geometric combination rule. The detailed scheme for parameterization will be discussed 
shortly. 
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Table 15. The Potential Parameters used for Na-Y, Ag-Y and H-USY 



Na-Y 






C-C 


H-H 


C-H 








A(kcal/mol) 


88371 


2861 


15901 








B(l/A) 


3.60 


3.74 


3.67 








C(A 6 *kcal/mol) 


583.13 


32.60 


137.88 










O-C 


O-H 


Na-C 


Na-H 






e(kcal/mol) 


0.25547 


0.16515 


0.04187 


0.03153 






a(A) 


2.9984 


2.5893 


3.3279 


2.7651 


Ag-Y 






O-C 


O-H 


Ag-C 


Ag-H 


(Case 1) 




e(kcal/mol) 


0.25547 


0.16515 


0.7823 


0.5895 






a(A) 


2.9984 


2.5893 


3.5773 


3.0145 


Ag-Y 


Uaz 




O-C 


O-H 


Ag-C 


Ag-H 


(Case II) 




e(kcal/mol) 


0.25547 


0.16515 


0.7601 


0.5728 






a(A) 


2.9984 


2.5893 


3.5773 


3.0145 


H-USY 






O-C 


O-H 










e(kcal/mol) 


0.1687 


0.1271 










G(A) 


2.9984 


2.5893 







1) Parameters for are the same for all sorbents. 

2) The numbers with underline were determined in this work. 
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Assignment of charges for adsorbate and adsorbent 

Assignment of charges also has a significant impact on simulation results 
when the Coulombic interactions are not negligible. In this work, the existing model for 
benzene in Cerius 2 was used. In this model, each carbon in benzene carries -0.093 point 
20 charge, while hydrogen carries +0.093. These values were nearly the average values 

from the various charge determination methods: MOP AC (0.1022), ZINDO (0.072) and 
charge equilibration method (Q eq ) by Rappe and Goddard (0.1 147) (Rappe and Goddard, 
1991). 

The point charges of zeolite were assigned either by (i) using the point 
25 charges in the literature or (ii) calculating by means of the charge equilibration method 
in the Cerius2 software. In this work, Na, Si/Al, O in Na-Y carried the charges of +0.78, 
+1.2125 and -0.72, respectively, since these values were successfully used in quantum 
mechanics simulations (Kitagawa et al., 1996). Here, no distinction was made between 
Si and Al atoms in terms of point charges. Charges of Ag-Y and all-silica Y-Zeolites 
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were calculated by the charge equilibration method: Ag(+0.45), Si/Al(+1. 16875), 0(- 
0.65), which will be referred to as Case I for Ag-Y, and Si(+1.2650), O(-0.6325) for H- 
USY. In order to investigate the effect of the assigned charges to the isotherm 
simulation results, the same charges that were used for Na-Y, i.e., Ag(+0.78), 
5 Si/Al(+1 .2125) and O(-0.72), were also used for Ag-Y. This case is referred to as Case 
II for Ag-Y. For all cases, regardless of the positions of cations and framework atoms 
(Si, Al, O), the same values of point charge were assigned to the same atoms throughout 
the simulation unit cell, even though the point charge might vary with the position of the 
atom. The Ag cation charges used for Monte Carlo simulation were also consistent with 
10 the values calculated by molecular orbital method at MP2/3-21G (Ag charge : +0.49) or 
MP2/LanL2DZ level (Ag charge : +0.74) using the zeolite cluster model (Chen and 
Yang, 1996). 

Results and Discussion 
15 Experimental 

Benzene/cyclohexane adsorption isotherms 

Pure component isotherms of benzene and cyclohexane were measured 
for H-USY(Si/Al=195), Na-Y(Si/Al=2.43), Ag-Y(Si/Al=2.43), and Pd-Y(Si/Al-2.43) at 
120 and 1 80 °C, shown in Figures 22-25. Isosteric heats of adsorption calculated using 
20 the Clausius-Clapeyron equation are listed in Table 16. 

Table 16. Comparison of Heat of Adsorption (kcal/mol) between Experimental and 
Simulated Results 







Na-Y 


Ag-Y 


H-USY 


Benzene 


Experimental 


17.7-18.2 


20.5-20.7 


6.2-7.2 




Monte Carlo Simulation 


16.4-18.5 


18.8-21.1 


8.2-9.2 




Molecular Orbital Calculation 


*** 


17.2 


*** 




Barthomeuf et al. 1} 


17.5-18.0 


*** 


*** 


Cyclohexane 


Experimental 


10.9-13.0 


12.1-12.3 


4.1-4.3 




Barthomeuf et al. 1} 


12.0-13.2 


*** 


*** 



1) Barthomeuf and Ha, 1973 



In the case of H-USY, the adsorbed amounts of cyclohexane were higher 
30 than that of benzene at > 0.04 atmosphere. This result qualitatively agreed with 
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calculations using the Horvath-Kawazoe theory (Horvath and Kawazoe, 1983; Cheng 
and Yang, 1994). The threshold pore-filling vapor pressures of cyclohexane and 
benzene at 120°C and 180°C were calculated using the spherical pore model of the H-K 
method. The threshold pressure is the pressure where the steep rise in the isotherm 
5 occurs. At both temperatures, it was found that the pore-filling pressure of cyclohexane 
was 20% lower than that of benzene. The reason for this difference was the larger 
polarizability (1 LOxlO" 24 cm 3 for cyclohexane and 10.3xl0" 24 cm 3 for benzene) and 
magnetic susceptibility (11 .3xl0~ 29 cm 3 for cyclohexane and 9. lxl 0" 29 cm 3 for benzene) 
of cyclohexane. However, the calculated isosteric heat of adsorption of benzene was 

10 slightly higher than that of cyclohexane. A possible reason for the larger heat of 

benzene adsorption was the existence of H + in the zeolite, although the number of proton 
was extremely small (less than one per unit cell). The interaction between protons in 
zeolite and 7i-electrons of benzene might be stronger than the interaction between 
oxygen in zeolite and benzene. Therefore, the amounts of benzene adsorption were 

15 slightly larger than that of cyclohexane at less than 0.04 atm, leading to its higher heat of 
adsorption. 

On all other sorbents investigated in this work, larger amounts of benzene 
were adsorbed at all pressure ranges. Heats of adsorption for benzene were 1 8 kcal/mol 
for Na-Y and 21 kcal/mol for Ag-Y, while heats of adsorption for cyclohexane were 10- 

20 12 kcal/mol for both Na-Y and Ag-Y. These results are clearly due to the strong 

interactions between the cations in zeolite and the benzene molecule. Consequently, 
Ag-Y and Na-Y could adsorb large amounts of benzene even at the pressure of less than 
0.01 atm. The heats of adsorption on Na-Y were comparable to previous work 
(Barthomeuf and Ha, 1973). Most importantly, Ag-Y showed the strongest interactions 

25 for benzene and the largest heats of adsorption (20.5 kcal/mol). For cyclohexane 
adsorption, the adsorption amounts and heat of adsorption were almost the same 
between Ag-Y and Na-Y, since Ti-complexation was absent for cyclohexane. A detailed 
discussion will be made in the section of molecular orbital calculation and Monte Carlo 
simulation results. These results indicated that improved purification performance of 

30 benzene from cyclohexane is possible using Ag-Y. The isotherms of both Ag-Y and 
Na-Y are reversible, although it took a little longer time to desorb benzene and 
cyclohexane. For example, 200 minutes are necessary to desorb benzene completely 
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after the adsorption at 0. 1 atmosphere. This slow desorption is clearly due to the strong 
interaction between benzene and Ag. However, the bonds are of the order of 20 
kcal/mol, which are weak enough so a substantial fraction of the adsorbed benzene is 
desorbed in a few minutes. For practical applications, thermal desorption may be a 
5 desirable option. 

Pd-Y gave lower adsorption amounts than both Ag-Y and Na-Y. This 
decrease of adsorption amounts by ion-exchange with Pd was due to degradation of the 
zeolite structure. In Figure 26, the pore size distributions for all sorbents calculated 
from N 2 isotherms at 77 K are shown. Pd-Y had a small pore volume with a wide pore 

10 size distribution, clearly indicating the collapse of the zeolite pore structure or pore 

blockage by Pd. The Pd ion-exchanged was performed at 80 °C, compared with the Ag 
ion-exchange at room temperature. Collapsing of zeolite structure tends to occur 
considerably more in acidic solution and at higher temperatures. The pore volume of 
Ag-Y also decreased by 23% after ion-exchange. However, this reduction of Ag-Y was 

15 mainly caused by 37% density increase of the zeolite lattice by replacing Na + with Ag + 
(27% reduction of pore volume per weight basis). The pore structure of Ag-Y was 
intact after ion-exchange. 

Another interesting phenomenon in Figure 26 is the apparent shift of the 
pore size distribution in H-USY toward a higher value, compared to other zeolites. This 

20 was likely due to the lower polarizability of oxygen in the framework of high-Si0 2 

zeolites (Pellenq and Nicholson, 1993), which will be discussed in detail in the section 
of Monte Carlo simulation results. Otherwise, the interactions between nitrogen and 
cations in zeolites would have led to a shift of the pore sizes toward smaller values for 
Na-Y, Ag-Y and Pd-Y. Since the pore sizes of Na-Y, Ag-Y and Pd-Y were the same as 

25 the actual value for Y zeolites (1 .3nm), the latter explanation is not plausible. 

Separation factors 

The pure component adsorption data were first fitted by both D-A 
(Dubinin- Astakhov) and L-F (Langmuir-Freundlich) isotherms. The separation factors 
30 of benzene over cyclohexane were then obtained using the equivalent multi-component 
isotherms. In the D-A equation of pure component gas, the volume adsorbed, V, at 
relative pressures P/P s is expressed as 
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V = V 0 exp 




(10) 



where 



RT 



(11) 



0Eo 



Doong and Yang (Doong and Yang, 1988) extended the D-A equation to mixed-gas 
adsorption in a simple way by using the concept of maximum available pore volume 
without any additional equations such as the Lewis relationship (Lewis et al., 1950). For 
binary mixtures, 



Since equation (12) and (13) are linear, V, and V 2 can be obtained in a straightforward 
fashion. After the calculation of V, and V 2 , volumetric adsorbed amounts were 
converted to molar adsorbed amounts and the separation factors were calculated. The 
conversion of molar adsorbed amount to volumetric adsorbed amount and vice versa 
was performed using the liquid density data of the adsorbates. The liquid density at 120 
and 180°C were obtained by the modified Rackett equation by Spencer and Adler 
(Spencer and Adler, 1978). The liquid densities used in this work were 0.769 g/cc at 
120°C and 0.689 g/cc at 180°C for benzene, and 0.679 g/cc at 120°C and 0.606 g/cc at 
180°C for cyclohexane. The adsorbed molar volumes were assumed to be the same in 
both single- and mixed-gas adsorption. 

For the Langmuir-Freundlich isotherm equation, pure component 
isotherm is written using three adjustable parameters: 




12) 



r 2 =(*"«- P.) exp 




(13) 




1 + BP 



(14) 
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The L-F equation can also be extended for n-component mixture (Yang, 1987): 



The fitting parameters for both isotherms are listed in Table 17 and the 
fitted curves are plotted in Figures 22 -25. The separation factors of benzene over 
cyclohexane were calculated using the equilibrium mole fractions in adsorbed and gas 
phases: 



X lY 

benzene I benzene 

CI ~ (16) 



X lY 

cyclohexan e / cyclohexan e 



The calculated separation factors for Ag-Y and Na-Y are compared in 

10 Table 18. Both D-A and L-F equations led to similar dependence of separation factor on 
pressure. The separation factor for Ag-Y increased with decreasing concentration of 
benzene, while that for Na-Y at 180° C decreased. Although the separation factors 
calculated by the two isotherms were not the same, the advantages of Ag-Y over Na-Y 
are clearly seen, especially at lower concentrations of benzene. This is certainly because 

15 Ag-Y can maintain its adsorption ability to benzene even at around lxlO" 5 atm. On the 
other hand, at 0.01 and 0.001 atm of benzene and at the lower temperature of 120°C, 
Na-Y showed equivalent separation factors to Ag-Y. 

For separations by adsorption and membranes, a separation factor of the 
order of 10 is considered high. For example, in air separation by pressure swing 

20 adsorption using zeolites and also by polymeric membranes, systems with separation 
factors (N 2 /0 2 for adsorption and 0 2 /N 2 for membranes) of near 10 are commercially 
used. For aromatics/aliphatics separations using liquid membranes (e.g., 
toluene/heptane and benzene/heptane separations), separation factors of the order of 10- 
100 yielded excellent separations (Li, 1968; Li, 1971). The separation factors shown in 

25 Table 18, in the order of 10-10,000, are excellent for separation/purification. 
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Table 17. Parameters for the D-A and L-F Equations 











D-A Equation 




L-F Equation 




Adsorbent 


Adsorbate 


Temp. (°C) 


v 0 


C 


n 


Ps 


q Q 


B 


n 








(cc/g) 


(-) 


(-) 


(atm) 


(mmol/g) 


(1/atm) 


(1/atm) 


Na-Y 


Benzene 


120 


0.295 


0.092 


3.76 


2.96 


3.15 


30.7 


2.53 




Benzene 


180 


0.326 




d 7^ 


10.07 


Z..OZ. 


449 


0.9 




Cyclohexane 


120 


0.267 


ft ]AQ 


S 01 


2.84 


? S8 


SI 2 


1.18 




Cyclohexane 


180 


0.264 


0 168 


6 44 


9.48 


2.30 


61.9 


0.77 


Ag-Y 


Benzene 


120 


0.244 


0.056 


1.57 


2.96 


3.07 


3.67 


7.46 




Benzene 


180 


0.248 


0 OR? 


1 OS 


10.07 


9 0^ 


1 S6 


1.85 




Cyclohexane 


120 


0.225 


0.133 


4.32 


2.84 


2.02 


48.4 


1.44 




Cyclohexane 


180 


0.237 


0 145 


4.96 


9.48 


1.80 


74.3 


0.91 


H-USY 


Benzene 


120 


0.237 


0.318 


2.84 


2.96 


3.40 


11.2 


0.62 




Benzene 


180 


0.234 


0.333 


2.17 


10.07 


3.40 


1.11 


0.74 




Cyclohexane 


120 


0.249 


0.265 


5.85 


2.84 


1.33 


2500000 


0.2 




Cyclohexane 


180 


0.15 


0.232 


4.93 


9.48 


0.73 


260 


0.39 


Pd-Y 


Benzene 


120 


0.173 


0.237 


2.15 


2.96 


0.883 


240 


0.7 




Cyclohexane 


120 


0.16 


0.317 


1.34 


2.84 


0.753 


9.22 


1.2 



20 Table 18. Separation Factors estimated by the D-A and L-F Equations 



30 



Pressure (atm) 




D-A Equation 


L-F Equation 


Benzene 


Cyclohexa 


Temp. (°C) 


Ag-Y 


Na-Y 


Ag-Y 


Na-Y 




ne 












0.01 




120 


107 


*** 


6 


12 


0.001 




120 


416 


441 


45 


48 


0.0001 




120 


2286 


1016 


335 


191 


0.00001 




120 


14361 


2553 


2462 


767 


0.01 




180 


39 


95 


19 


5.4 


0.001 




180 


120 


71 


57 


4.2 


0.0001 




180 


456 


26 


163 


3.3 


0.00001 




180 


1767 


0.5 


470 


2.6 
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Diffusion time constants 

In Table 19, diffusion time constants calculated from 0 - 0.3 fractional 
uptake were summarized. Diffusion time constants of benzene and cyclohexane in Na- 
Y and Ag-Y were of the order of 10" 4 1/s. These values were almost the same as those of 
5 benzene and cyclohexane into AgN0 3 monolayer dispersed Si0 2 sorbent (Takahashi et 
al.). However, these values were an order of magnitude lower than those of propylene 
and propane in AgN0 3 monolayer-dispersed Si0 2 sorbent (Rege, 1998) and at least the 
same as those of ethylene and ethane into Ag ion-exchanged resins (Wu et al., 1997). 
The diffusion rates of benzene in these zeolites are high for applications (Yang, 1987). 

10 



Table 19. Diffusion Time Constants of Benzene and Cyclohexane 



Adsorbents 


Temp. 


Benzene 


(Pressure Change) 


Cyclohexane 


(Pressure 




(°C) 


(1/s) 


(atm) 


(1/s) 


Change) 












(atm) 


Na-Y 


120 


2.1xl0" 4 


(5.4xl0- 3 -->2,6xl0' 2 ) 


9.8X10- 4 


(l.lxlO" 2 - 


Ag-Y 


120 


4.2x1 0" 4 


(5.4xl0' 3 -->2.6xl0- 2 ) 


l.lxlO- 4 


>2.0xl0' 2 ) 


AgNO 3 /SiO 2 (0.33g/g) '> 


120 


1.3xl0- 4 


(2.0xl0' 3 «>3.7xl0- 2 ) 


2.9xl0 4 


(l.lxlO 2 -- 












>2.0xlO' 2 ) 












(2.0xl0- 3 - 












>3.7xl0 2 ) 


Adsorbents 


Temp. 


Propylene 


(Pressure Change) 


Propane 


(Pressure 




(°C) 


(1/s) 


(atm) 


(1/s) 


Change) 












(atm) 


AgNO 3 /SiO 2 (0.33g/g) 2 > 


70 


3.5xl0" 3 


(0~>1) 


*** 


*** 


Ag-Amberlyst-35 3) 


70 


1.5xl0" 4 


*** 


L4xl O" 4 


*** 



1 ) Takahashi et al. , 2000 
20 2)Regeetal, 1998 
3) Wu etal., 1997 

Molecular Orbital Calculation 
Adsorption Bond Energy and Molecular Geometries 
25 Adsorption energy and structural geometry were calculated using 

molecular orbital theory. The calculated heat of adsorption using B3LYP/3-21G level 
by equation (1) is also compared with experimental and Monte Carlo simulation results 
in Table 16. The theoretical value, 17.2 kcal/mol, was in fairly good agreement with the 



61 

experimental value. The comparison of bond length by benzene adsorption revealed that 
the Ag-O length increased from 3. 023 A to 3.246A. The C-C bond length in benzene 
also increased from 1.385A to 1.385-1. 391 A. This trend was consistent with that for 
benzene interaction with AgCl (Takahashi et al.). 

5 

Natural Bond Orbital (NBO) Results 

The nature of benzene- Ag zeolite interaction can be better understood by 
analyzing the NBO results. The main feature of the interaction can be seen from the 
population changes in the vacant outer-shell s orbital of the silver and those in the d 

10 orbitals of silver upon adsorption. In Table 20, NBO analysis results for benzene on Ag 
zeolite are compared with other combinations of adsorbates and adsorbents. According 
to the electron population changes by benzene adsorption on Ag zeolite, the trend is in 
agreement with previously reported trend of 7i-complexation (Cheng, 1996; Huang et al., 
1999a; 1999b; 1999c ; Takahashi et al.) i.e., donation of electron charges from the p- 

15 orbitals of benzene to the vacant s-orbital of the silver (a donation) and, simultaneously, 
back-donation of electron charges from the d-orbitals of silver to 7t* orbital of benzene 
(d-7i* back-donation). The contribution of d-7i* back-donation is larger than a donation, 
resulting in the decrease of net electron charge in silver. The net electron charge of 
carbon in benzene increased instead. 

20 The differences of electron population changes between benzene and 

ethylene adsorption deserve a discussion. From Table 20, it is seen that the contribution 
of donation is predominant in ethylene adsorption, while d-7c* back-donation is slightly 
larger in benzene adsorption. Also, the extent of a donation and d-7i* back-donation by 
benzene adsorption was considerably smaller compared with ethylene adsorption. 

25 Therefore, 7t-complexation with benzene was weaker than that with ethylene. Weak n- 
complexation of benzene can also be understood by the difference in bond lengths. The 
bond length between carbon in ethylene and silver in AgCl was 2.66 A (Huang et al., 
1999b), while that between carbon in benzene and silver in Ag zeolite or AgCl was 
2.71-3.7 A or 3.22 A (Takahashi et al.), respectively. This longer bond length was 

30 caused by the less overlap of orbitals from the weaker 7t-complexation. 
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Table 20. Summary of the NBO Analysis of the n-complexation: Electron Population Changes of 
Each Orbital After Adsorption 









Silver 






Carbon 






Method/Basis Set 


5S 


S4d 


Net Change 


2s 


S2p 


Net Change 


Benzene on Ag Zeolite 


B3LYP/3-21G 


0.00541 


-0.0135 


-0.00805 


0.00554 


0.00709 


0.01263 


Benzene on AgCl 0 


B3LYP/LanL2DZ 


0.0033 


-0.0071 


-0.0038 


*** 


*** 


*** 


C 2 H 4 on Ag Zeolite 2) 


MP2/3-21G 


0.0596 


-0.0177 


0.0419 


0.0203 


0.0193 


0.0396 


C 2 H 4 on AgCl 2) 


MP2/3-21G 


0.1204 


-0.0470 


0.0734 


0.0226 


0.0142 


0.0368 


C 2 H 4 onAgCl 3) 


B3LYP/LanL2DZ 


0.061 


-0.055 


0.006 


*** 


*** 


*** 



1) Takahashi et al., 2000 

2) Chen and Yang, 1996 

3) Huang etal., 1999 



Monte Carlo Simulation 

Monte Carlo simulation of adsorption isotherms provides great insight 
into the adsorption mechanism via model construction, charge assignment and forcefield 
determination. In this section, benzene adsorption on Na-Y, H-USY and Ag-Y is 
simulated and the mechanism of adsorption on each zeolite is discussed. 

Na-Y zeolite 

GCMC simulation results of benzene isotherms at 120 and 180°C using 
the potential parameters given in Table 15 are shown in Figure 27, where they are also 
compared with experimental data. The simulation results were in excellent agreement 
with experimental data. The calculated heats of adsorption were 16.4 - 18.5 kcal/mol, 
which were in agreement with that from experiments. It was hence confirmed that the 
energy expressions and parameters were suitable for simulation of benzene isotherms on 
Na-Y. 

H-USY zeolite 

First, simulation of benzene isotherms on H-USY was attempted using 
the same potential parameters as on Na-Y. However, the simulated isotherms were 
considerably higher than the experimental data. This disagreement was caused by 
incorrect well depth parameters, e 0 _ c and e Q _ H , for van der Waals interaction. Pellenq 
and Nicholson (Pellenq and Nicholson, 1993) reported the values of the polarizabilities 
of framework atoms in silicalite and aluminosilicates determined from Auger electron 
spectroscopy data. It was found that the polarizability for oxygen in zeolite framework 
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is sensitive to the ratio of silicon to aluminum and also the nature of the extraframework 
cations. Watanabe et al. (Watanabe et al., 1995) used different potential values for 
oxygen in silicalite (0.202 kcal/mol), denoted by Oz(Si-O-Si), and that in 
aluminosilicates (0.334 kcal/mol), denoted by Oz(Si-O-Al). 
5 In this work, an adjustment of e 0z(Si . 0 . si) . c and e 0z(Si . 0 . Si) . H was made to 

simulate the experimental Henry's law constants at 120 and 180°C. Here since low 
pressure data depend only on the benzene-zeolite interactions so the benzene-benzene 
interactions were neglected. Using the well depth parameter for e Na _ c = 0.04187 
kcal/mol and e Na . H = 0.03153 kcal/mol for Na-Y (in Table 15) and assuming e Na = 0.041 

10 kcal/mol as used by Watanabe et al. (Watanabe et al., 1995) for Na-A and Na-X, e c = 
0.0428 kcal/mol and e H = 0.0243 kcal/mol were obtained from the geometric 
combination rule, equation (9). Therefore, £ 0 z(si-o-si)-c e oz(Si-o-so-H f° r H-USY could 
be expressed as (0.0428 x e 0z(Si _ 0 _ Si) ) 0 - 5 and as (0.0243 x e 0 z(Si-o-si))° 5 > respectively. By 
fitting with e 0z (si-o-si)> oz(Si-o-so = 0-665 kcal/mol was obtained to fit the experimental 

15 Henry's law constant of 0.519 molecules/(cell.kPa) at 120°C and 0.137 

molecules/(cell.kPa) at 180°C. The ratio of e 0z (si-o-si) over 8 oz(si-o-/vi) was 0.44 (calculated 
from c 0z(si . 0 . A1) _ H ) - 0.60 (calculated from e 0z (si-o-Ai)-c )> which was in agreement with the 
ratio of 0.60 obtained by Watanabe et al. (Watanabe et aL, 1995). 

Using the modified parameters for H-USY (in Table 15), the benzene 

20 isotherms were calculated and compared with experimental isotherms in Figure 28. 
Again, the isotherms and the heat of adsorption (8.2 - 9.2 kcal/mol) obtained from 
simulation were in excellent agreement with the experimental data. 
Ag-Y zeolite 

The potential parameters for Ag-Y are listed in Table 15. These 
25 parameters were modified to get the low pressure points of the isotherm at 1 80°C (1 .30 
mmol/g at L6 x 10" 4 atmosphere), although adjustment of the parameters to fit the 
experimental Henry's law constant was preferable. The reason not to use the 
experimental Henry's constant was that the small pressure values ( 8.7 x 10' 6 
atmosphere) at the lowest adsorbed amount in the isotherms might lead to the large error 
30 of potential parameters. The procedure of adjustment was basically the same as in the 
case for H-USY except for changing the e Ag value instead of the e 0z(Si . 0 . Si) value. e c = 
0.0428 kcal/mol and e H = 0.0243 kcal/mol were obtained from the geometric 
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combination rule, equation (9). Then, e Ag . c and e Ag _ H for Ag-Y could be estimated from 
(0.0428 x 6 Ag ) 0 5 and (0.0243 x e Ag )° \ respectively. In Case I for Ag-Y, e Ag became 14.3 
kcal/mol, while e Ag = 13.5 kcal/mol was obtained for Case II. The larger well-depth 
parameter was needed for Case I to compensate for the weaker Coulombic interaction 
5 due to the smaller Ag cation charge. In both cases, the e Ag values used were more than 2 
orders of magnitude greater than e Na , clearly indicating the strong interaction from k- 
complexation. 

Figure 29 shows the simulation results of benzene on Ag-Y for Case I 
and Case n. The differences between Case I and Case II were fairly small. In both cases, 
10 the agreement with experimental isotherms was not as good as that for Na-Y and H- 

USY. However, they certainly followed the basic trend of the experimental isotherms. 
The calculated heat of adsorption was 18.8-21.1 kcal/mol, which was higher than Na-Y 
and also in agreement with the experimental value. 

In all cases, not only Ag-Y but also Na-Y and H-USY, the simulated 
1 5 isotherms tended to be higher than the experimental isotherms in the lower pressure 
range with few exceptions, while they were lower in the higher pressure range. The 
higher simulation results at low pressures might be caused by the larger values for the 
well-depth parameter between oxygen (zeolite) - carbon (benzene) interaction and 
oxygen (zeolite) - hydrogen (benzene) interaction. The other possible cause was the 
20 imperfection of the zeolite crystals. Since zeolites contain defects while in Monte Carlo 
simulation perfect zeolite lattice was assumed. In the high pressure range, larger 
experimental values were possibly caused by adsorption between zeolite crystals. 

Adsorption sites for benzene on Na-Y and Ag-Y were investigated using 
the Mass Cloud Plots function in Cerius2 (not shown here). It was confirmed that 
25 benzene was adsorbed preferentially near the SIL cation sites and the center of the 12- 
ring window, which was consistent with the neutron diffraction data for Na-Y obtained 
by Fitch (Fitch et al., 1986). For Ag-Y, benzene was adsorbed slightly further away 
from the SII cation sites than that in Na-Y, due to the fact that the van der Waals radius 
of Ag + is 30% larger than that of Na + . 



30 
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Notation 

A = repulsive parameter for Buckingham potential 

B = repulsive parameter for Buckingham potential or Langmuir constant 

C = dispersive parameter for Buckingham potential or constant for D-A equation 

D = diffusivity 

E = energy 

n = constant 

P = pressure 

q = point charge or molar adsorbed amount 

R = van der Waals radii or gas constant 

r = distance or radius 

T = temperature 

U = potential energy 

V = volumetric adsorbed amount 

X = equilibrium mole fraction in adsorbed phase 

Y = equilibrium mole fraction in gas phase 
Greek letters 

a = polarizability or separation factor 
(3 = adjustable parameter 

e = well depth parameter for Lennard- Jones potential 
a = distance parameter for Lennard- Jones potential 
Subscript 

i = atom site or component 
j = atom site or component 
s = saturation 
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Figs. 30-32 are graphs depicting further breakthrough curves. 

5 EXPERIMENT E 

Adsorbent Preparation. The adsorbents used in this study were copper 

cation forms of Y Zeolite and also type PCB Activated Carbon (PCB-AC). The former 

were prepared initially using conventional liquid phase ion exchange techniques. The 

starting material, sodium type (Na) Y Zeolite (Strem Chemicals, Si/Al = 2.43), was 

10 used as received and in powder form. The activated carbon samples were obtained from 
Calgon Corporation also in powder form and used without further purification. 

Cu(I)-Y (or reduced Cu(IT)-Y) was prepared by first ion exchanging Na-Y 
with a Cu(N0 3 ) 2 aqueous solution (0.5 M) for 48 hours followed by reduction of Cu 2+ to 
Cu + . The amount of copper in the ion exchange solution was equivalent to a 5-fold 

1 5 cation exchange capacity. The adsorbent was recovered by filtration and washed with 
copious amounts of deionized water, followed by drying at 100°C for at least 24 hours. 
Activation was performed at 450°C in pure helium to promote auto-reduction of Cu 2+ 
species to Cu + , which is desired for p-complexation (Cu(I)-Y). Auto-reduction of cupric 
ions to cuprous ions in synthetic zeolites has been reported by several groups (Larsen, 

20 S.C.; Aylor, A.; Bell, A.T.; Reimer, J.A. Electron Paramagnetic Resonance Studies of 
Copper Ion-Exchanged ZSM-5. J. Phys. Chem. (1994), 98, 1 1533; Parrillo, D.J.; 
Dolenec, D.; Gorte, R. J.; McCabe, R.W. Adsorption Studies on CuZSM-5: 
Characterization of the Unique Properties of Ion-Exchanged Cu. J. Catal. (1993), 142, 
708 and Iwamoto, M.; Yahiro, H.; Tanda, K.; Mizuno, N.; Mine, Y.; Kagawa, S. 

25 Removal of Nitrogen Monoxide through a Novel Catalytic Process: 1. Decomposition 
on Excessively Copper Ion-Exchanged ZSM-5 Zeolites. J. Phys. Chem. (1991), 95, 
3727) and a more detailed discussion about the Cu-Y auto-reduction process is 
discussed earlier. 

30 Reagents and Standards: Gasoline and diesel samples were obtained from British- 
Petroleum (BP) gas stations located in the Michigan area. Gasoline was unleaded 
regular type, while diesel is type #2 according to BP's specifications. The average total 
sulfur concentration for these fuels was reported to be 335 and 430 ppmw, respectively 
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(data available from BP). Thiophene, benzothiophenes (BT), and dibenzothiophene 
(DBT) standards were purchased from Sigma- Aldrich. 

Elemental Analysis, Cu(II)-Y and Na-Y adsorbents were characterized using neutron 
activation analysis (NAA) in the research nuclear reactor of the Phoenix memorial 
Laboratory at the University of Michigan. The data was obtained from a 1 -minute core- 
face irradiation delivered via pneumatic tube to a location with an average thermal flux 
of 2.13 x 10 12 neutrons/cm 2 /s. Following irradiation, two separate counts were made, 
one after a 13 -minute decay (for Al and Cu) and a second count after a 1 hour and 56 
minutes decay (for Na and K); both were for 500 seconds. 

The concentrations of Al, Na, and K were determined based on comparison 
with three replicates of the standard reference material NIST1633A (coal fly ash); the 
determination of Cu content was evaluated relative to high-purity copper wire. Data 
reductions for NIST1633A were based on NIST certified values. 

Fixed-Bed Adsorption/Breakthrough Experiments. All adsorption/breakthrough 
experiments were performed in vertical custom made quartz adsorbers equipped with a 
supporting glass frit. The setup consisted of a low-flow liquid pump equipped with a 
ceramic piston and cylinder liner, Kynar compression fittings, two Pyrex feed tanks, and 
a heating element. Initially, the adsorbents were loaded inside the adsorber (between 1 
or 2 grams), and heated in situ (450°C) while flowing helium. Temperatures were 
maintained at steady values using a PK) temperature controller. The gases used for 
activation were pretreated inline prior to contacting the sorbent using a 3 A-type zeolite. 
The latter allows removal of trace water from the gases, which may adsorb in the 
adsorbents that are being tested. After activation treatment, the adsorbent under study 
was allowed to cool down to room temperature under helium and then tapped to ensure 
proper packing. Next, a sulfur- free hydrocarbon was allowed to flow downwards 
through the sorbent at a rate of 0.5 cm 3 /min. After wetting the adsorbent for about 10 
minutes, the feed was changed to either commercial grade gasoline or diesel at a 0.5 
cmVmin rate. Samples were collected at regular intervals until saturation was achieved, 
which depended on the adsorption dynamics and amount of adsorbent. 
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Gas Chromatography Analysis. All the samples collected during the breakthrough 
experiments were analyzed using a Shimadzu GC-17A v3 unit equipped with an EC-5 
capillary column (L = 30 m; ID = 0.32 mm) and a flame photometric detector (FPD). 
The column temperature program was set to increase from 50°C to a set value at a 
5 5°C/min rate. For gasoline and diesel sulfur analysis, the column final temperature 

values were 250°C and 330°C, respectively. In addition, a split mode injection was used 
for all samples at a 100:1 ratio. About 4mL of sample volume was injected for each GC- 
FPD run. The injection port temperature was set to either 250°C or 330°C, this 
depending on the nature of the samples under analysis. 

10 Peak identification information for sulfur compounds present in gasoline 

and diesel was gathered after using standards and by retention time comparison with 
data available in the literature. 5 A Ma, X. et al., Catal Today. (2002), 77, 107; Ma, X. et 
al.; Am. Chem. Soc. Div. Pet. Chem. Prepr. (2002), 47, 48; Gates, B.C.; Tops0e, H. 
Reactivities in Deep Catalytic Hydrodesulfurization Challenges, Opportunities and the 

1 5 Importance of 4-Methyldibenzothiophene and 4,6-Dimethyldibenzothiophene. 

Polyhedron (1997), 16, 3213; Bacaud, R.; Cebolla, V.L.; Membrado, L.; Matt, M.; 
Pessayre, S.; Galvez, E.M. Evolution of Sulfur Compounds and Hydrocarbons Classes 
in Diesel Fuels during Hydrodesulfurization. Combined Use of Thin-Layer 
Chromatography and GC-Selective Chemilumiscence Detection. Ind. Eng. Chem. Res. 

20 (2002), 41, 6005 and Ma, X.; Sakanishi, K.; Isoda, T.; Mochida, I. Determination of 
Sulfur Compounds in Non-Polar Fraction of Vacuum Gas Oil. Fuel. (1997), 76, 329. 
For standards, thiophene, BT and DBT solutions were diluted in sulfur free ^-octane to a 
known concentration and then injected for retention time determination. The total sulfur 
concentration during breakthrough adsorption experiments for either gasoline or diesel 

25 was estimated after integrating the entire GC chromatogram region. It was assumed the 
total concentration is proportional to the collective areas under the peaks and that it 
follows a linear fashion. This is believed to be acceptable for conditions where 
complete sample elution after injection is achieved and after correction for noise data. 
The total elution times for gasoline and diesel for the GC conditions described above 

30 were 30 and 50 minutes, respectively. A similar procedure was used to calculate 

individual sulfur components concentration during adsorption experiments. For this 
purpose, only the specific compound unique peak area was used. Figures 33 and 34 
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show GC-FPD chromatograms for as-received gasoline and diesel fuels. The results are 
similar to those found in literature for similar analysis conditions and set-up. 



FTIR Spectroscopy Analysis. Fuels were analyzed for aromatic/aliphatic C-H bond 
5 stretching using FTIR spectroscopy on a Nicolet Impact 400 FTIR spectrometer 

equipped with a TGS detector. The samples were loaded into a liquid IR cell fitted with 
ZnSe windows prior each analysis following standard procedures. The spectra were 
then taken at room temperature using 100 scans per run and a resolution of 4 cm* 1 . 
Background spectra were collected using the liquid IR cell without any sample inside. 

10 

Results and Discussion 

Adsorbent Characterization. Characterization for the adsorbents used in this study 
was achieved by NAA. All the zeolites were in hydrated conditions before testing. As 
seen in Table 21, ion-exchanging Na-Y with Cu 2+ results in an incomplete ion exchange. 

15 If it is assumed that one Cu 2+ cation compensates for two aluminum tetrahedra, then in 
the prsent case, the ion exchange resulted in only 70% substitution of the original 
sodium ions (i.e., 2Cu/Al = 0.72). The remaining sodium species were then 
compensating for the other aluminum tetrahedra; in other words, the (2Cu + Na)/Al ratio 
should be unity, which is observed here. This scheme is the simplest one used to 

20 describe such behavior. More elaborated schemes can be found elsewhere. 2530 Parrillo, 
et al., J. Catal 1993, 142, 708 and Baes, C.F.; Mesmer, R.E. The Hydrolysis of Cations; 
Wiley; New York, 1976. 

Adsorbent Activation, Copper Auto-reduction and Migration. Since Y zeolite is 
25 known to be highly hydrophilic (uptake > 20 wt% water at ambient conditions), all the 
gases used for activation of the adsorbents were pretreated with 3A-Type zeolites prior 
to entering the fixed-bed unit. For Cu(II)-Y, the activation was performed at 450°C in 
helium to promote auto-reduction of Cu 2+ species to Cu + , which is desired for n- 
complexation (Cu(I)-Y). Larsen et al. reported that about 40% of the cupric ions in 
30 CuZSM-5 are reduced under helium at 410°C. Larsen, S.C. et al., J, Phys. Chem. 1994, 
98, 1 1533. Also, EPR studies done by Takahashi et al. showed that 50% of Cu 2+ in 
Cu(II)-Y zeolite is reduced under vacuum or helium at 450°C, (Takahashi, A.; Yang, 
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R.T.; Munson, C.L.; Chinn, D. Cu(I)-Y-Zeolite as a Superior Adsorbent for 
Diene/Olefin Separation., Langmuir (2001), 17, 8405) which is in good agreement with 
the findings of Larsen et al. It should be mentioned that after 18 hours of auto- 
reduction, the color of Cu(I)-Y was pale green, compared to a bluish green typically 
5 observed in Cu(II)-Y. This was sure evidence of auto-reduction, as Cu + should result in 
a white color (like in CuCl). In the present invention, successful liquid phase thiophene 
adsorption experiments as described hereinabove with Cu(I)-Y zeolites provide further 
evidence of reduction of Cu 2+ . 

The adsorption capacity of auto-reduced Cu-Y observed for simple 

10 liquid hydrocarbon mixtures is believed to be due not only to the reduction of the copper 
ions, but also because of the final position of the cations after activation and/or during 
adsorption. Cu + cations occupy exposed sites in the Y zeolite framework, such as sites 
II and III according to the nomenclature of Smith, (Smith, J.V. Faujasite-Type 
Structures. Aluminosilicate Framework. Positions of Cations and Molecules. 

15 Nomenclature. Adv. Chem. Ser. (1971), 101, 171) in order to interact with the sulfur 

molecules. Recent studies by Fowkes et al. have shown that, upon reduction of Cu(II)- 
Y, there was a redistribution of cation positions and most of the reduced species (Cu + ) 
occupied sites I* and II. Fowkes, A. J.; Ibberson, R.M.; Rossensky, M.J. Structural 
Characterization of the Redox Behavior in Copper-Exchanged Sodium Zeolite Y by 

20 High Resolution Powder Neutron Diffraction. Chem. Mater. (2002), 14, 590. Lamberti 
et al. showed similar results for both reduced Cu(H)-Y and Cu(I)-Y prepared by ion 
exchange and gas phase reaction with CuCl, respectively. Lamberti, C; Spoto, G.; 
Scarano, D.; Paze, C; Salvalaggio, M.; Bordiga, S.; Zecchina, A.; Turnes Palomino, G.; 
D' Acapito, F. Cu l -Y and Cu n -Y Zeolites: a XANES, EXAFS and visible-NIT Study. 

25 Chem. Phys. Lett. (1997), 269, 500. In addition, some Cu + ions could be induced to 

migrate to more exposed sites under the presence of guess molecules. Turnes-Palomino 
et al. found, based on IR data, that this happens when CO molecules are adsorbed in 
Cu(I)-Y. Tumes Palomino, G.; Bordiga, S.; Zecchina, A.; Marra, G.L.; Lamberti, C 
XRD, XAS, and IR Characterization of Copper-Exchanged Y Zeolite. J. Phys. Chem. B 

30 (2000), 104, 8641 . Their zeolite was also prepared by vapor-phase exchange of H-Y 
with CuCl. Thus, based on Turnes-Palomino 's findings, there could also be some 
synergy that contributes for the observed adsorption capacity of this adsorbent. 
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Fixed-Bed Adsorption Tests. After activating the adsorbents in-situ, a solution of 
either commercial regular gasoline or diesel was passed through the fixed-bed and the 
downstream sulfur concentration monitored as a function of time. Figure 35 shows a 
breakthrough/adsorption profile obtained during desulfurization of gasoline with a 
Cu(I)-Y adsorbent. Adsorption capacities at both breakthrough and saturation were 
obtained after integrating for the area above the breakthrough curves. For gasoline 
treatment with Cu(I)-Y the loadings at breakthrough and saturation were 0.14 and 0.39 
mmol/g, respectively (Table 22). In the present invention, breakthrough loadings were 
observed of 1.82 and 0.22 mmol/g after thiophene (2,000 ppmw) removal from w-octane 
and benzene/w-octane (20 wt% C 6 H 6 ), respectively. These results indicate that the 
adsorbent performance for sulfur removal is greatly affected by the presence of 
aromatics, which should occur during gasoline treatment. Moreover, since other 
aromatic species besides benzene can still interact with 7t-complexation adsorbents, it is 
important to consider the whole composition spectrum for gasoline. Table 23 shows 
typical aromatic composition for both the commercial gasoline and diesel. Gasoline has 
a considerable amount of aromatic molecules larger than benzene (because of additional 
functional groups), and they should compete for void space during adsorption. It should 
be mentioned that King et al. found a similar result for vapor phase sulfur removal from 
toluene mixtures with ZSM-5 zeolites. King, D.L.; Faz, C; Flynn, T. Desulfurization of 
Gasoline Feedstocks for Application in Fuel Reforming. SAE Paper 2000-01-0002; 
Society of Automotive Engineers: Detroit, MI, 2000. They concluded the aromatics 
entities were responsible for a diminishing in adsorbent uptake capacity. 

Tournier at al. studied the adsorption of xylene isomers in liquid phase on a K,Ba-Y 
zeolite. Tournier, H.; Barreau, A.; Tavitan, B.; LeRoux, D.; Sulzer, C; Beaumont, V. 
Two Experimental Methods to Study Adsorption Equilibria of Xylene Isomers in the 
Liquid Phase on a Y Zeolite. Micropor. Mesopor, Mat. 2000, 39, 537. For a mixture of 
m-xylene and p-xylene, K,Ba-Y adsorbs 0.13 mmol/g of the latter for an equilibrium 
concentration (liquid phase) is 2.4 mole%. This provides additional evidence that 
aromatics will compete for adsorption sites, even within themselves. However, even in 
the presence of compounds such as xylenes, Cu(I)-Y still offers considerable selectivity 
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towards the sulfur in gasoline and this has to be due to 7i-complexation. One gram of 
adsorbent is capable of producing 14 cm 3 of sulfur free gasoline (< 1 ppmw S), 
removing not only thiophenes and benzothiophenes (BTs), but also their substituted 
counterparts. 

5 

Figure 36 depicts the evolution of each sulfur species downstream Cu(I)-Y as 
determined by GC-FPD analyses. Breakthrough curves for selected gasoline sulfur 
compounds (Figure 37) were obtained also from the GC-FPD data. The chromatograms 
show the adsorbent has selectivity towards substituted thiophenes and heavier molecules 

10 over non-substituted molecule. This can be observed in the relative intensities of the 
peaks in the 4 minutes (elution time) region (please refer to Figure 33 for peaks 
reference). It should be mentioned most substituted thiophenes and BTs remain highly 
unreacted during conventional HDS processes, which is clearly not the case when using 
7r-complexation zeolites. 

1 5 One alternative to increase both breakthrough and saturation capacities is 

to provide a guard bed to take care of the aromatics during the 7r-complexation process 
with Cu(I)-Y. PCB type activated carbon (AC) was chosen to accomplish this. PCB- 
AC is used commercially for both liquid- and vapor-phase applications including 
recovery of alcohols, hydrocarbons and aromatics. The present inventors performed 

20 vapor-phase single component equilibrium adsorption experiments for thiophene and 

diesel in PCB- AC. More benzene than thiophene was adsorbed even at pressures below 
1 x 10 2 atm, which indicates polarizability plays an important role in the performance of 
this adsorbent. A similar behavior should be observed also in liquid-phase adsorption. 
Choma et al. studied how to predict solute adsorption behavior form liquid-phase 

25 adsorption on carbons based on the corresponding gas/solid adsorption parameter, and 
vice- versa. Choma, J. Burakiewicz,-Mortka, W.; Jaroniec, M.; Gilpin, R.K. Studies of 
the Structural Heterogeneity of Microporous Carbons Using Liquid/Solid Adsorption 
Isotherms. Langmuir 2000, 39, 537. They demonstrated, for both types of interfaces, 
that similar information about the structural heterogeneity of the adsorbent can be 

30 gathered. It should be mentioned that their results were obtained for benzene as solute. 

For the experiments involving activated carbon, the guard bed was 
loaded on top of the Cu(I)-Y bed after making sure the latter was packed uniformly. The 
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adsorbent activation procedure was identical to the one described in the experimental 
section. Figure 35 shows the resulting sulfur breakthrough curve for gasoline treatment 
after adding an activated carbon (type PCB) guard bed on top of Cu(I)-Y. The guard 
bed accounted for 15 wt% of the total adsorbent bed. For this experiment, the new 
5 adsorbent combination proves to enhance the sulfur removal process, indicating this the 
activated carbon probably offers less mass transfer resistance towards aromatics 
allowing for an apparent sulfur concentration wave front to move ahead faster and reach 
Cu(I)-Y first. During the momentary absence of aromatics, the sulfur compounds 
should take advantage of the 7i-complexation process since the competitive adsorption 

10 process observed for "unprotected" Cu(I)-Y is delayed. Breakthrough and saturation 

capacities (from Figure 35) were 0.18 and 0.50 mmol/g, respectively (Table 22). This is 
a 28% increase in capacity over the one offered by Cu(I)-Y without a guard bed. 

It is to be understood that the guard bed may be located in any suitable 
position. In a non-limitative embodiment, the guard bed is located adjacent the inlet of 

1 5 the main bed. It is to be further understood that the guard bed may be formed from any 
suitable sorbent, including but not limited to activated carbon, activated alumina, silica 
gel, zeolites, clays or pillared clays, diatomaceous earth, any of the porous sorbents, 
and/or mixtures thereof. 

GC-FPD analysis (Figure 38) shows that thiophene molecules 

20 breakthrough faster, possibly indicating that competitive adsorption now lies (at least for 
some time) within the sulfur compounds (substitutes and non-substituted) themselves. 
This can be seen also in Figure 39, where thiophene concentration has reached almost 
60% of the corresponding saturation value after processing about 30 cm 3 of gasoline per 
gram of adsorbent. The methyl-substituted BTs have reached only 10% on their 

25 respective concentration saturation values at similar conditions. 

The same adsorbent combination (AC/Cu(I)-Y) was also tested for diesel 
desulfurization. Figures 40, 41, and 42 show total sulfur breakthrough curves, GC-FPD 
chromatograms, and selected sulfur containing molecules breakthrough curves, 
respectively. For total sulfur removal from diesel, the breakthrough and saturation 

30 loadings are 0.32 and 0.59 mmol/g, respectively. This translates to a saturation capacity 
of about 1 .85 wt% sulfur. The GC-FPD data and individual component breakthrough 
curves show the adsorbent combination has some selectivity towards specific molecules. 
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Figure 42 shows 4,6-dimethyldibenzothiophene molecules leaving the bed early, but not 
in considerable fashion. This behavior could be due to pore clogging with some 
substituted dibenzothiophenes. However, this is not expected since Y zeolite is known 
to adsorb substituted aromatics molecules in considerable quantities. For example, at 
5 P/P 0 = 0.5, X zeolite (FAU, Si/Al = 1.33) adsorbs 43.3 benzene molecules and about 37 
toluene molecules per unit cell, respectively. Na-X (Si/Al = 1.25) also adsorbs larger 
molecules such as (C 4 H 9 ) 3 N (-16 molecules per unit cell), which has a kinetic diameter 
of a = 8. 1 A. Breck, D.W. Zeolite Molecular Sieves, Wiley: New York, 1973, reprinted 
by Krieger, Malabar, Fla., 1984. Thus, pore clogging in Cu(I)-Y sure occur only for 

10 heavily substituted dibenzothiophenes and in the presence of framework faults, specially 
those affecting pore window locations. 

When comparing treatment results for both gasoline and diesel and the 
same adsorbent combination (i.e. Cu(I)-Y with a guard bed) it is believed that the fuel 
composition (other than the sulfur content) plays an important role in the adsorbent 

1 5 performance. Since, as discussed above, the amount of aromatics is believed to play an 
important role in diminishing the sulfur removal capacity (i.e. aromatics will also do k- 
complexation with copper cations), differences in aromatic content between different 
fuels should be taken into consideration. Although Table 23 shows some aromatics 
such as benzene as not as abundant in diesel than in gasoline, an estimate of total 

20 aromatic content is in place. 

C-H bond stretching for both aromatics and aliphatics is sometimes 
straightforward to identify from FTIR spectroscopy. C-H stretching for alkyl, alkenyl, 
and aromatics shows intensities in the 2853 - 2962, 3010 - 3095, and -3030 cm" 1 
regions, respectively. Graham-Solomons, T. W. Fundamentals of Organic Chemistry, 

25 2nd ed.; Wiley: New York, 1986. JR intensities are usually strong for the first and third 
regions, while for the second case these are of a moderate type. Figure 43 presents FTIR 
spectra for the gasoline and diesel used during this study. Stahl and Tilotta got similar 
results for gasoline and a lighter fuel. Stahl, D.C.; Tilotta, D.C. Partition Infrared 
Method for Total Gasoline Range Organics in Water Based on Solid Phase 

30 Microextraction. Environ. Sci. Technol (1999), 33, 814. From a qualitative perspective, 
it is clear the C-H bond stretching relative intensities are stronger near the aromatic 
range for gasoline, while diesel shows a stronger signal in the aliphatic region. An 
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estimate of the aromatic/aliphatic composition was obtained after integrating the spectra 
for the regions just mentioned. The total aromatic content for gasoline and diesel was 
estimated to be 33.7 and 24.7%, respectively. Thus, since linear molecules (such as 
alkanes) are not expected to compete during adsorption, less aromatics should favor 
5 sulfur removal with Cu(I)-Y or AC/Cu(I)-Y. 

The present inventors were able to completely regenerate Cu(I)-Y 
zeolites after saturating them with thiophene molecules removed from liquid 
hydrocarbons. They found that using air at 350°C provides a suitable way to completely 
remove the sulfur containing molecules completely. After this treatment, the adsorbent 
10 needs to be auto-reduced in an inert atmosphere to recover its activated form. It is 

believed the same regeneration scheme may be used after treating gasoline or diesel, but 
the adsorbent adsorption capacity is enough to consider the use of disposable units for 
on-board applications. 

15 Conclusion 

Results have demonstrated that copper (auto-reduced) Y-type zeolites are excellent 
adsorbents for removal of all sulfur compounds form commercial liquid fuels, based on 
fixed-bed adsorption experiments. Adsorbents were capable of reducing sulfur content 
to values less than about 1 ppmw sulfur for long periods of time. When used with a 
20 guard bed, Cu(I)-Y provides the best adsorption capacity both at breakthrough point and 
at saturation. It is believed the guard bed delays the aromatics concentration wave front 
and this allows for exclusive adsorption of sulfur entities with Cu(I)-Y zeolites. 

Table 21. Composition data for Cu(II)-Y and Na-Y zeolites obtained from neutron activation analysis (NAA). 
25 Molar Ratios 

Adsordent Na/Al Cu/Al Si/Al 

Na-Y 0.94 - 2.43 

Cu(II-Y) 0.28 0.36 2.43 
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Table 22 .Breakthrough and saturation loadings for gasoline or diesel total sulfur in 
Cu(I)-Y zeolites. 



Fuel Type 


Average Sulfur Content 


Adsorbent 


Breakthrough Loading 


Saturation Loading 




(wt%) 




(mmol/g) 


(mmol/g) 


Gasoline 


335 


Cu(I)-Y 


0.14 


0.39 






AC/Cu(I)-Y" 


/LIS 


0.50 


Diesel 


430 


AC/Cu(I)-Y" 


CU4 


0.59 



Adsorbent bed made of two layers: top is PCB type activated carbon (AC) and bottom is Cu(I)-Y 



10 zeolite. AC accounts for 15 wt% of the bed weight. 



Table 23. Typical aromatic composition of some commercial fuels. Average values 
1 5 are shown in parentheses. 



Typical Composition (vol%) 



Fuel 


Total 


Benzene 


Toluene 


Xylene 


Ethylbenz 


Cumene 


Type 


Aroma tics 








ene 




Gasoline 


12.8-56.7 
(30.7) 


0.4- 
3.48 
(1.16) 


1.5-16.8 
(7.0) 


0.8-8.4 
(4.0) 


0.36-2.63 
(1.36) 


0.04-0.17 
(0.08) 


Diesel 


15.7-47.2 
(31.8) 


0 - 0.02 
(0.01) 


0-0.18 
(0.05) 


0.04-0.71 
(0.15) 


0.02-0.7 
(0.04) 


0.01-0.06 
(0.02) 



Note: Data available form British Petroleum (BP), 2001 TRW Summer Motor Gasoline Survey Report, 
and AAM Diesel Fuel Surveys. 



25 EXPERIMENT F 

After activating the adsorbents in-situ, gasoline was passed through the 
fixed-bed, and the sulfur concentration in the effluent was monitored as a function of 
time. Fig. 44 shows the breakthrough curve obtained during desulfurization of gasoline 
with a Cu(I)-Y adsorbent. Adsorption capacities at both breakthrough (i.e., the break 

30 point, or the first point with detectable sulfur) and saturation were obtained upon 

integration of the area above the breakthrough curve. For gasoline treatment with Cu(I)- 
Y, the loadings at breakthrough and saturation were 0.14 and 0.39 mmol/g, respectively 
(Table 24). From Fig. 44, one gram of Cu(I)Y was capable of producing 14 cm 3 of 
sulfur free gasoline (< 0.28 ppmw S), removing not only thiophenes and 

35 benzothiophenes (BTs), but also their substituted counterparts. Fig. 44 also shows the 
resulting sulfur breakthrough curve for gasoline after adding an activated carbon (type 
PCB) guard bed on top of Cu(I)-Y. The guard bed accounted for 15 wt% of the total 
adsorbent bed. For this run, the new adsorbent combination proved to significantly 
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enhance the sulfur removal capacity. The breakthrough and saturation capacities (from 
Fig. 44) were 0.18 and 0.50 mmol/g, respectively (Table 24). 

Fig. 45 depicts the evolution of each sulfur species in the adsorber 
effluent as determined by GC-FPD analyses. The chromatograms show the adsorbent 
5 has selectivity towards substituted thiophenes and heavier molecules over non- 
substituted molecules (i.e., thiophene). This can be observed from the relative 
intensities of the peaks in the 4-minute (elution time) region. 

Table 24. Adsorption Capacity for Desulfurization of Commercial Gasoline 

10 

mmol/g wt% S 

Cu(I)-Y Breakthrough 0.21 0.67 

Capacity 

Saturation 0.59 1.85 

Capacity 

Activated Breakthrough 0.27 0.87 

15 Carbon Capacity 

(15%)/Cu(I)-Y 

Saturation 0.74 2.31 

Capacity 

Note: A total initial sulfur concentration of 500 ppmw was assumed for the above quantities 

20 EXAMPLES OF SOME SUITABLE ADSORBENTS 

A brief description of some non-limitative examples of adsorbents 
which may successfully be used in the present invention follows. Detailed 
descriptions may be found in U.S. Patent No. 6,423,881, and in U.S. Patent No. 
6,215,037, each of which patents is incorporated herein by reference in its entirety. 

25 The adsorbent comprises a carrier having a surface area, the carrier 

having dispersed thereon a monolayer of a d-block transition metal compound, 
eg., a silver compound. The silver compound releasably retains the thiophenes; 
and the carrier comprises a plurality of pores having a pore size greater than the 
effective molecular diameter of the thiophenes. 

30 It is to be understood that any suitable carrier may be used. 

However, in a preferred embodiment, the carrier has a BET surface area greater 
than about 50 square meters per gram and up to about 2,000 square meters per 
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gram, and comprises a plurality of pores having a pore size greater than about 3 
angstroms and up to about 10 microns. In a more preferred embodiment, the 
carrier is a high surface area support selected from the group consisting of 
refractory inorganic oxide, molecular sieve, activated carbon, and mixtures 
5 thereof. Still more preferred, the carrier is a refractory inorganic oxide selected 
from the group consisting of pillared clay, alumina and silica. 

It is also to be understood that any suitable silver compound may 
be used. However, in a preferred embodiment, the silver compound is a silver (I) 
halide. In a more preferred embodiment, the silver compound is a silver salt, and 

10 the salt is selected from the group consisting of acetate, benzoate, bromate, 
chlorate, perchlorate, chlorite, citrate, fluoride, nitrate, nitrite, sulfate, and 
mixtures thereof. 

In one exemplary embodiment of this embodiment of the present 
invention, the silver compound is silver nitrate (AgN0 3 ) and the carrier is silica 

15 (Si0 2 ). 

The method of the present invention may further comprise the step 
of changing at least one of the pressure and temperature to thereby release the 
thiophenes-rich component from the adsorbent. It is to be understood that the 
pressures and temperatures used may be within a suitable range. However, in the 

20 preferred embodiment, the selected pressure of preferential adsorption is a first 
pressure, and the pressure of release is a second pressure less than the first 
pressure. In a more preferred embodiment, the first pressure is in a range of about 
1 atmosphere to about 35 atmospheres, and the second pressure is in a range of 
about 0.01 atm to about 5 atm. 

25 In the preferred embodiment, the selected temperature of 

preferential adsorption is a first temperature, and the temperature of release is a 
second temperature greater than the first temperature. In a more preferred 
embodiment, the first temperature is in a range of about 0°C to about 50°C, and 
the second temperature is in a range of about 70°C to about 200°C. 

30 Without being bound to any theory, it is believed that the retaining 

of the thiophenes is accomplished by formation of 7c-complexation bonds between 
the silver compound and the thiophenes. Separation by 7c-complexation is a 
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subgroup of chemical complexation where the mixture is contacted with a second 
phase, which contains a complexing agent. The advantage of chemical 
complexation is that the bonds formed are stronger than those by van der Waals 
forces alone, so it is possible to achieve high selectivity and high capacity for the 
5 component to be bound. At the same time, the bonds are still weak enough to be 
broken by using simple engineering operations such as raising the temperature or 
decreasing the pressure. 

The 7i-complexation generally pertains to the main group (or d- 
block) transition metals, that is, from Sc to Cu, Y to Ag, and La to Au in the 

10 periodic table. These metals or their ions can form the normal a bond to carbon 
and, in addition, the unique characteristics of the d orbitals in these metals or ions 
can form bonds with the unsaturated hydrocarbons (olefins) in a nonclassic 
manner. This type of bonding is broadly referred to as n -complexation, and has 
been considered for gaseous hydrocarbon separation and purification using 

1 5 cumbersome liquid solutions. 

A further example of an adsorbent comprises an ion-exchanged 
zeolite selected from the group consisting of zeolite X, zeolite Y, zeolite LSX, and 
mixtures thereof, the zeolite having exchangeable cationic sites, and at least some 
of the sites have a d-block transition metal cation, eg., silver cation or copper 

20 cation, present. In the more preferred embodiment, the cationic sites of the ion- 
exchanged zeolite contain silver cation. In a further preferred embodiment, 
essentially all cationic sites of the ion-exchanged zeolite contain the silver cation. 

Without being bound to any theory, it is believed that the 
preferential adsorption occurs by 7i-complexation. 

25 In an embodiment of the present invention, gasoline has been 

desulfurized at about 24 cc/g sorbent. Using a standard cartridge size that contains 
5 kg of an embodiment of the present inventive sorbent, approximately 120 liters 
of gasoline may be desulfurized. For a fuel cell automobile, this advantageously 
translates into a driving range of about 7200 miles (or about 120 hours). The 

30 sorbent of embodiments of the present invention may also be regenerated for 
reuse, although this is not necessary, as the inventive sorbent is relatively 
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inexpensive. For a conventional internal combustion engine automobile, the 
resultant driving range is approximately 3600 miles. 

While preferred embodiments of the invention have been described 
in detail, it will be apparent to those skilled in the art that the disclosed 
5 embodiments may be modified. Therefore, the foregoing description is to be 

considered exemplary rather than limiting, and the true scope of the invention is 
that defined in the following claims. 



